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Résumé 


Les semi-conducteurs magnétiques dilués (DMS) sont des matériaux 
promoteurs pour les futures applications spintroniques. Nous les avons 
investigués, dans le cadre de cette thèse, en dopant le semi-conducteur binaire 
MgeSe par des éléments de transition tels que le manganèse (Mn) et le 
chrome (Cr). Leurs principales propriétés physiques ont été obtenues par des 
calculs de premier principe en se basant sur la théorie de la fonctionnelle de 
la densité, en spin polarisé (SP-DFT) et sur la méthode des ondes planes 
augmentées et linearisées avec un potentiel complet (FP-LAPW), qui est 


implémentée dans le code WIEN2K. 


Les composés qui ont fait objet de cette thèse admettent la formule générale 
MeixIM,Se (TM= Mn, Cr et x = 0.25, 0.50 et 0.75). Nous avons calculé leurs 
propriétés structurales, électroniques et magnétiques. Les calculs des 
propriétés thermodynamiques n’ont concerné que les composés dopés au 


chrome (Cr). 


Le calcul des propriétés structurales des alliages (Ma:,Mn,Se et Ma ,Cr,Se) a 
été effectué en utilisant l’approximation WC-GGA ; les résultats obtenus ont 
servi au calcul des propriétés thermodynamiques des composés Mo1.Cr,Se. 
Quant aux propriétés électroniques et magnétiques, c’est l’approximation 


TB-mBJ combinée avec Papproximation PBE-GGA qui à été utilisée. 


Les résultats ont montré que les composés Mo:.Mn,Se sont des semi- 
conducteurs ferromagnétiques pour les deux directions de spins (up et down) 


tant dis que les composés Moi,Cr,Se sont des demi-métaux. 


Les mots clé : Premier-principe, MoSe, Semi-conducteurs magnétiques 


dilués, Spintronique. 


Abstract 


Diluted magnetic semiconductors (DMS) are promoter materials for future 
spintronic applications. We have investigated them, as part of this thesis, by 
doping the binary semiconductor MgeSe with transition elements such as 
manganese (Mn) and chromium (Cr). Their main physical properties were 
obtained by first-principle calculations based on the spin-polarization density 
functional theory (SP-DFT) under the framework of the full potential 
linearized augmented plane waver (FP-LAPW) which is implemented in the 
WIEN2K code. 


The compounds which were the subject of this thesis have the general 
formula Me: ,TM,Se (TM= Mn, Cr and x= 0.25, 0.50 and 0.75). We calculated 
their structural, electronic, and magnetic properties. The calculation of 
thermodynamic properties concerned only the compounds doped with 


chromium (Cr). 


Calculation of structural properties of alloys Moe:.Mn,Se and Mer, Cr,Se) was 
performed using the WC-GGA approximation; the results obtained were 
used to calculate the thermodynamic properties of the Ma,CrSe 
compounds. As for the electronic and magnetic properties, the TB-mB] 


approximation, combined with the PBE-GGA approximation, was used. 


The results showed that the Mgi,Mn,Se compounds are ferromagnetic 
semiconductors for the two types of spins (up and down) as long as the 


compounds Mo:.Cr,Se are half-metals. 


Keywords: First-principle, MoSe, Diluted magnetic semiconductors, 


Spintronic. 
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INTRODUCTION 
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Les propriétés électroniques et magnétiques des matériaux sont exploitées 
dans les technologies de l'information modernes. Comme aux premiers temps 
de l'informatique, le traitement de l'information dans les ordinateurs et toutes 
les machines numériques, tablettes, téléphones portables, etc .…., est effectué 
dans des circuits à l'état solide où les états logiques "1" et "0" sont associés au 
potentiel électrique, et sont représentés par des transistors ou positions de 
mémoire (bits) qui prennent l'une ou l'autre valeur selon le courant électrique 


reçu. 


Contrairement à l’électronique classique, le but de l'électronique de spin 
(spintronique) [1] est de créer de nouveaux composants électroniques qui 
exploitent à la fois la charge de l’électron ainsi que son moment de spin. De 


tels systèmes, permettent de stocker un grand volume d'informations. 


L’une des avancées majeures dans le domaine de la spintronique est la 
magnétorésistance géante (GMR) [2], découverte en 1988 par A. Fert et P. 
Grünberg et couronnée par le prix Nobel de la physique en 2007. Ce 
phénomène a été observé pour la première fois dans un système multicouche, 
formé d’un empilement de couches alternées de Fer (Fe) et de Chrome (Cr) 
selon la configuration (Fe/Cr/Fer/Substrat) [3]. La conductance d’un tel 
système dépend fortement du champ magnétique appliqué, qui peut 
commuter le couplage entre les couches formant l’empilement. Ce 
phénomène (GMR) s'applique dans le dispositif à vannes de spin [4] et il est 


utilisé dans les têtes de lecture de disques durs dans les ordinateurs [5]. 


Dans ce contexte, la réalisation de matériaux avec des propriétés à la fois semi- 
conductrices et magnétiques à température ambiante (ou à haute température) 
reste toujoufrs d'actualité et constitue un défi pour la communauté des 


chercheurs dans le domaine de la physique des matériaux. A la fin des années 
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1970, Jaczynski et al (1978) [6], et Gaj et al (1978) [7] ont mis une stratégie 
pour créer de tels systèmes en introduisant des moments magnétiques locaux 
dans les matrices de semi-conducteurs. Cette stratégie a donné naissance à 


une nouvelle classe de matériaux, appelés semi-conducteurs magnétiques 


dilués, en abrégé (DMS) [8]. 


Depuis, les recherches en spintronique ont été principalement axées sur les 
semi-conducteurs magnétiques dilués à trois dimensions (les semi- 
conducteurs classiques II-VI, III-V ...etc.) mais depuis la découverte du 
graphène en 2004 [9], elles ont été étendues aux systèmes à deux dimensions 
comme : le Nitrure de bore, les dichalcogénures de métaux de transition 
(TMD), nitrure de bore hexagonal (h-BN), Cristaux à boucle monoatomique 


émergents ‘Xenes’...etc.) 


Des recherches approfondies continuent dans ce domaine, elles ont comme 
objectif la magnétisation de matériaux semi-conducteurs par des procédés tels 
que le dopage par substitution [10-12], par l'adsorption d'atomes en surface 


[13] et par création des défauts (des lacunes) dans le système [14]. 


Dans ce travail de thèse, nous étudions leffet du dopage par substitution sur 
les propriétés du semi-conducteur MgSe qui constitue la matrice hôte. Ce 
composé est parmi les éléments les plus importants de la famille des semi- 
conducteurs de type II-VI, en raison des avantages qu’il présente, comme 
nous le ferons dans cette thèse. Les dopants choisis sont des éléments 
magnétiques : le chrome (Cr) et le manganèse (Mn) qui possèdent tous les 


deux des électrons de type ‘d’. 


Pour se faire, nous calculons les différentes propriétés de nos composés en se 
basant sur la théorie de la fonctionnelle de la densité (DFT), qui est une 
méthode de calcul ab-initio et elle est de nos jours la plus appropriée pour 


prédire le plus fidèlement possible les propriétés des solides. 
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D'une manière générale, les outils de la simulation numérique ne cessent de 
se développer pour une meilleure prise en charge de la caractérisation des 
matériaux se présentant sous formes multiples. Ces outils ont l’avantage d’être 
indépendants des mesures expérimentales et ne prennent en compte que les 
caractéristiques de base (masse, nombre de charge, structure atomique, etc.) 


des éléments considérés. 


Les résultats délivrés par les calculs ab-initio sont dans la plupart du temps en 
parfait accord avec les résultats obtenus par les études expérimentales. L’un 
des exemples est la demi-métallicité, connue dans le domaine de la 
spintronique. Cette caractéristique a été découverte pour la première fois avec 
des calculs théoriques (ab-initio) [15], et elle a été confirmée plus tard 


expérimentalement [16]. 
Le travail que nous présentons dans cette thèse comprend quatre chapitres : 


Dans le premier chapitre, nous présenterons un aperçu sur le rôle du spin 
dans l'électronique et quelque généralité sur les semi-conducteurs 


magnétiques dilués. 


Le second chapitre est consacré au formalisme de la théorie de la 
fonctionnelle de la densité (DFT). Les notions de base et les différentes 
approximations utilisées dans le cadre de cette théorie sont présentées. Nous 
avons également présenté dans les détails la méthode des ondes planes 
augmentées linéarisées (FP-LAPW) ainsi que le code de calcul WIEN2K que 


nous avons utilisé. 


Dans le troisième chapitre, nous avons discuté les résultats des propriétés 
structurales, électroniques et magnétiques du composé MoSe dopé au 


manganèse (Mn) dans la structure zinc-blende (ZB). 
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Dans le chapitre quatre, nous avons étudié les mêmes propriétés mentionnées 
précédemment du composé MeSe dopé au Chrome (Cr) dans la structure 


rock-salt (RS), avec en plus les propriétés thermodynamiques. 


Enfin, une conclusion englobant les résultats les plus importants est donnée 
avec des perspectives ouvrant la voix sur l’investigation d’autres propriétés 
comme complément aux résultats réalisés dans le cadre de cette thèse de 


doctorat. 
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L.1 Introduction 


Dans cette partie de thèse, nous allons dans un premier temps nous intéresser 
à l’électronique de spin, appelée aussi « spintronique » qui est une nouvelle 
branche de la physique, située à l’interface entre l'électronique et le 
magnétisme. Le concept général de cette branche est d’utiliser le spin de 
l’électron, comme second degré de liberté, en plus de sa charge, pour traiter 


le transport dans les solides. 


Nous allons tout d’abord présenter en revue l’histoire du magnétisme dans 
l’électronique et notamment la découverte de la magnétorésistance géante et 
de la magnétorésistance tunnel, ainsi que les domaines d’applications de 
lélectronique de spin et les bouleversements technologiques importants 
qu’elle à créés dans le domaine de la microélectronique. Nous donnerons 
ensuite un aperçu sur les nombreuses directions de la spintronique et sur la 
nouvelle classe de matériaux qu’elle regroupe à l’image des semi-conducteurs 
magnétiques dilués (DMS) qui ont fait l’objet d’une grande partie des travaux 


de recherche dans le cadre de cette thèse de doctorat. 


Enfin, nous présenterons les interactions magnétiques observées dans les 


DMS ainsi que l'origine du ferromagnétisme dans ce type de matériaux. 
I.2 Historique de la Spintronique 


La microélectronique a réellement commencé en 1947 avec la découverte de 
l'effet transistor par les trois physiciens américains ]. Bardeen, W. Brattain et 
W. Shockley [1] qui ont obtenu, suite à cette invention, le prix Nobel en 1956. 
Cependant, le premier transistor à effet de champ n’a fonctionné qu’en 1959 


(par Atalla) [2] ; cette invention à booster l’électronique qui, jusque-là, s’est 
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basée sur les porteurs de charges (électron ou trous) pour coder et traiter 


l'information. 


Par ailleurs, les progrès technologiques accomplis en matière de composants 
microélectroniques ont connu une croissance exponentielle. Gordon Moore, 
un des fondateurs d’intel [3], a parlé du principe du doublement, tous les deux 
ans, du nombre de transistors présents sur une puce. Cependant, cette course 
à la miniaturisation a atteint aujourd’hui ses limites, et il devient de plus en 
plus difficile de continuer à améliorer de manière significative les 


performances des composants électroniques. 


Dans lélectronique classique, la fonction de stockage et le traitement des 
informations restent toujours physiquement et technologiquement séparés ; 
l'information est véhiculée par les courants de charge dans les canaux de 
conduction, alors que le stockage des données est réalisé par des composants 
magnétiques, cela nous amène directement à une nouvelle branche appelée 
électronique de spin ou encore la spintronique [4]. C’est une nouvelle branche 
de la physique, en plein expansion ; elle est basée sur la combinaison de la 
charge de l’électron et de son spin, chose qui a été négligée par l’électronique 


classique. 


L'utilisation du spin dans le stockage de l'information conduit à des 


avantages tel que : 
(1) La non volatilité de l’information stockée quand on éteint l'ordinateur. 


(2) L’obtention des vitesses de traitement des bits d'information beaucoup 


plus rapides. 


(3) La baisse significative de la consommation de électricité (un 
processeur standard utilise énormément de courant et doit-être 


refroidit efficacement). 


(4) Le stockage d’une forte densité d'informations. 
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La découverte en 1921 du spin grâce à l'expérience de Stern et Gerlach a mis 
en évidence le moment magnétique propre appelé moment de spin. Il prend 
deux valeurs opposées : +1/2 (désigné par up, vers le haut) et -1/2 (désigné 
par down, vers le bas). Il est responsable des propriétés magnétiques des 
matériaux solides. Dans un matériau non magnétique (métallique, isolant ou 
semi-conducteur), le spin ne joue aucun rôle par ce que le nombre d’électrons 
de spin up et de spin down est le même, ce qui donne un moment magnétique 
total nul. Par contre, dans les matériaux magnétiques, en plus de la charge 
électrique, le spin des électrons peut être exploité comme facteur important 
des porteurs électriques où le nombre des électrons de spin up et de spin 
down est différent ; en effet, c’est cette différence qui donne lieu à Papparition 
d’un moment magnétique macroscopique non nul. L’asymétrie de spin des 


porteurs du courant est le premier principe de base de la spintronique. 


La spintronique [4] désigne l’ensemble des phénomènes physiques qui 
s'intéressent aux propriétés des courants de spin dans les matériaux, ainsi qu’à 
leur control en la présence ou en l’absence d’un champ magnétique externe et 


ceci dans le but d’améliorer les capacités de stockage de l'information. 


Un de ces phénomènes physiques est l’anomalie observée par Nevill Mott [5] 
en 1936 sur la résistivité des métaux ferromagnétiques au voisinage de la 
température de Curie. Cette anomalie a été expliquée par le rôle primordial du 


spin dans le transport. 


Cette idée à été reprise par Albert Fert et Peter Grünberg, dans les années 70, 
quand ils ont montré expérimentalement l’influence du spin sur la mobilité 
des électrons dans des métaux magnétiques et ont abouti sur la découverte de 
la magnétorésistance géante (GMR) [6,7] puis la magnétorésistance à effet 
tunnel (TMR) [8,9]. Ces découvertes ont été à l’origine de beaucoup 


d'applications dans les domaines des mémoires non volatiles MRAM), des 
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composants pour la logique, des composants radiofréquences pour les 


télécommunications, et des capteurs de champs magnétiques hypersensibles. 


L’électronique de spin a été à l’origine d’un nouveau concept mettant en 
valeur le spin de l’électron à côté de sa charge électrique. Elle à conduit à la 
découverte d’une nouvelle classe de matériaux regroupant à la fois les 
propriétés semi-conductrices et magnétiques. Ces matériaux sont appelés 
semi-conducteurs magnétiques dilués (DMS : dluted magnetics semi-conducators). 
La spintronique continue à se développer dans de très nombreuses directions 
pour répondre à des besoins technologiques dans divers domaines 


d'applications. 
L.3 La magnétorésistance géante GMR 


La découverte de la magnétorésistance géante (GMR : Geant Magnetic 
Resistance) a eu lieu en 1988 [6,7,10] par deux équipes de chercheurs : l’équipe 
d'Albert Fert et celle de Peter Grünberg ; cette découverte a donné naissance 


à la spintronique. 


Ce n’est qu'après le développement des techniques de dépôt sous ultravide, 
exemple de l’épitaxie par jets moléculaires (EJM) au milieu des années 80, qu’il 
est devenu possible de déposer des multicouches ultrafines, notamment les 


multicouches magnétiques. 


Ainsi, le phénomène de la magnétorésistance géante a été observé, la premier 
fois, dans des multicouches composées de couches ultraminces de fer (Fe) 
ferromagnétiques, séparées par de très fines couches de chrome (Cr) non 
magnétiques Fe/Cr/Fe [6,7] et dont l’épaisseur ne dépasse pas quelques 
couches atomiques. À cette échelle de la matière (échelle nanométrique), c’est 
la physique quantique qui dicte ses lois et on se tourne vers l'exploitation du 


spin au lieu de la charge comme c’est le cas de l'électronique. 
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La figure [1, est un schéma illustrant le mécanisme général de la 
magnétorésistance géante dans des structures multicouches ferromagnétiques 
alternant deux couches F1 et F2 séparées par une couche non magnétique 


d'épaisseur inférieure au libre parcours moyen des électrons. 


À gauche, dans la configuration magnétique antiparallèle (AP) figure I.1 (a), les 
électrons circulent plus facilement quand leur aimantation pointe dans le 
même sens que l’aimantation des couches. Dans ces conditions, les électrons 
sont tous bloqués au début ou à la fin et la résistance de la bicouche est donc 
grande figure [.1(b) Quand on applique un champ magnétique externe, on 
inverse l’aimantation d’une couche et on tombe dans ces conditions sur la 
configuration magnétique parallèle (P) figure I[.1(a). Dans ce cas-là, la moitié 


des électrons ne sont plus bloqués et la résistance et bien plus faible. 


4 
R/R(H=0) 
(AP) (P) 


Spin FM NM FM Spin FM NM FM 


, 
{Fe 30 À/ Cr18A),, 


Æ 80 % 


e e 
(Fe 30 À / Cr 12 A). 


(Fe30 À / Cr9 À), 


. à : i i à : 
Re R R 40 30 -20 “10 O0 10 20 30 4 


[a] Champ magnétique {kG) 


Figure I.1 : (a) Schéma du mécanisme de la GMR et (b) Premières 
observations de la magnétorésistance géante à Orsay sur des multicouches 


Fe/Cr(001) [6]. 


En résumé, la GMR représente la différence relative en pourcentage entre les 
valeurs mesurées des résistances selon les configurations magnétiques 
parallèle et antiparallèle, elle fut nommée « géante » car elle peut atteindre des 


valeurs importantes. Elle à la forme suivante : 
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Rap-Rp (1) 
Rp 


GMR = 


Où RP et RypSsont les résistances équivalentes du système dans les 


configurations parallèle et antiparallèle respectivement : 


- 2R.r (12) 
PU R+r 

R+r 12 

Rap —= 2 4 ) 


Cette découverte, couronnée par un prix Nobel de physique en octobre 2007 
[11,12], a attiré l’attention, aussi bien sur les possibilités d'applications que sur 
l'intérêt fondamental d’expériences de transport électronique exploitant le 


spin. 


Ensuite, beaucoup d’efforts ont été consentis en matière de recherches 
expérimentales et théoriques pour l'exploitation de la GMR, dans le domaine 
de la spintronique et la plus connue de ses applications concerne les têtes de 


lecture des disques dur d’ordinateurs [13,14]. 
I.4 La magnétorésistance tunnel TMR 


La magnétorésistance à effet tunnel est aussi dans des multicouches de type 
F1/NM/PF2 ressemblant à celle de la GMR sauf que cette fois-ci on a deux 
couches ferromagnétiques séparées par une fine couche d’un matériau isolant 
(d'épaisseur de l’ordre du nanomètre) constituant une barrière d’énergie pour 
les électrons de conduction figure [.2. Dans ce cas, le passage du courant entre 


les deux électrodes ferromagnétiques se fait par effet tunnel. 


Couche (a) : couche ferromagnétique « libre » 
Couche (b) : barrière tunnel isolante 
Couche (c) : couche ferromagnétique « piégée » 


Figure I.2 : Empilement simple d’une jonction tunnel magnétique. 
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Dans le modèle de la TMR, la probabilité de transmission des électrons à 
travers la barrière isolante dépend de l'angle formé par le moment de spin des 
électrons et l’aimantation de l’électrode qui reçoit les électrons ayant traversé 
la barrière. De ce fait, quand on applique entre les deux électrodes de la 
jonction une tension électrique, le courant à travers celle-ci va dépendre de 
l’angle entre les aimantations des deux électrodes. Cet effet avait été observé 
pour la première fois en 1975, par le Français, Michel Jullière, à basse 
température dans des jonctions magnétiques semi-conductrices à base de 
germanium [15]. D’après la figure IL.3, la magnétorésistance tunnel est définie 
comme la variation relative de la résistance du système parallèle (PR) et du 


système antiparallèle (AP) des aimantations. 


Figure L.3 : (a) Jonction tunnel magnétique dans la configuration parallèle, 
les états disponibles au niveau de Fermi dans la première électrode trouvent 
des états disponibles dans la seconde électrode (b) Jonction tunnel 
magnétique dans la configuration antiparallèle, les électrons présents à 


l'énergie de Fermi ne trouvent pas d’états disponibles pour les accueillir. 


"Dans la configuration parallèle : les électrons up et down au niveau de 
Fermi vont trouver la même densité d’états de l’autre côté de la barrière 


tunnel. Le courant tunnel est donc important. 


"Dans la configuration antiparallèle : les électrons up auront une densité 
plus importante alors que les électrons down auront une densité plus 


faible. Il en résulte un courant tunnel plus faible. 
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On peut écrire alors la TMR : 


Rap — Rp _ 2P;P} (L.4) 


TMR = = 
RP 1—P,P, 


P, et P, sont les polarisations des électrodes 1 et 2. Cette formule s’applique si 
les deux électrodes ferromagnétiques possèdent la même polarisation de 


spin P. 


Mais l'intérêt pour les jonctions tunnel n’a vraiment pas connu d’essor 
qu'après les travaux expérimentaux de l’équipes de Jagadeesh Moodera en 
1995, effectués sur des jonctions amorphes à base d’alumine [8] réalisées à 
température ambiante et qui présentent des amplitudes de magnétorésistance 


trois fois supérieures à la GMR, de l’ordre de 30 % à 70%. 


Le phénomène de la jonction tunnel suscite un intérêt considérable du point 
de vue fondamental et applications notamment pour le stockage de 


Pinformation, les composants logiques ou les capteurs magnétiques. 
L.5 Développements applicatifs de électronique de spin 


Dès la découverte de la magnétorésistance géante (GMR) et la 
magnétorésistance à effet tunnel (TMR), le transfert de l’électronique de spin 
vers les applications a été très rapide. Dans cette partie, on donnera un aperçu 
suf les développements applicatifs de la spintronique dans plusieurs 


domaines. 
L.5.1 Têtes de lecteur de disques durs d’ordinateurs 


Un disque dure figure [.4 est composé d’un disque métallique qui peut tourner 
jusque à 15 mille Tours par minute et d’une tête de lecture/écriture fixée au 
cout d’un bras. Le disque métallique recouvert d’une très fine couche 
magnétique est un alliage de cobalt dans lequel chaque grain représente un 


minuscule aimant. 
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En mode écriture, le courant qui parcourt la bobine crée un champ 
magnétique qui polarise tout une zone de grain. Chaque zone matérialise un 
bit! d’information 0 ou 1 selon le sens du courant qui traduit le langage binaire 


de l'ordinateur. 


Secteur 


Bloc de 
4 secteurs 


Tête de lecture (GMR)/(TMR) 


li d'écriture inductive 


Domaines aimantés Média d'enregistrement 


Figure I.4 : (a) Schéma de principe d’un disque dur et (b) d’une tête de 


lecture/écriture magnéto résistive. 


En 1997, les disques durs ont été équipés d’une double tête (figure L.5), la tête 
d'écriture inductive classique et la tête de lecture contenant un capteur qui 
intègre soit un élément GMR soit un élément TMR, donc basé sur le 
phénomène de magnétorésistance. Ces têtes sont indépendantes dans leur 
fonctionnement, mais collées lune à l’autre et fabriquées simultanément sur 
le même substrat qui forme le patin de vol. La tête GMR/TMR permet de 
détecter un fin champ magnétique généré par une centaine de grains de cobalt 


seulement. 


TUn bit c’est une unité élémentaire d’information électrique pouvant prendre deux valeurs : 0 (il n'y a 
pas de courant), ou 1 (il y a du courant) 
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ner 
isolant 


En BU. 


Nano-aimant 


Couche TMR 


Ecran inférieur 


Figure I.5 : (a) Schéma et photo d’une tête GMR et (b) d’une tête TMR. Le 


courant de lecture (1) est indiqué par des flèches rouges. 


Le premier disque dur commercialisé par IBM en 1956, est constitué de 50 
disques de 61 cm de diamètre, il peut stocker 5 millions d’octets (=5 Mo ; un 
octet =8bits), pour un poids dépassant la tonne. Aujourd'hui, plus de 300.000 
milliards d'octets (300.000 Go) de données peuvent être accumulées sur cette 
même surface et la quasi-totalité des têtes de lecture pour disques durs utilise 
la GMR/TMR, permettant des densités de plus de 20 milliards de bits 
(gigabits) par centimètre carré de disque : l'équivalent d'environ 2500 romans 
par centimètre carré. Plus de 600 millions d'ordinateurs équipés de telles têtes 


de lecture sont produits chaque année. 
I.5.2 Les capteurs de champ magnétique 


Les capteurs magnétiques constituent un domaine où l’électronique de spin a 
permis d'apporter une nouvelle richesse, en amenant la possibilité d’une 
grande sensibilité sur une petite échelle spatiale basé sur le phénomène de la 
GMR ou TMR. Ils sont implantés dans un très grand nombre d’objets 
technologiques qui nous entourent : dans un réfrigérateur, une machine à 


laver et même certaines chaussures de sport. 
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Les capteurs de champ magnétique ont de nombreuses applications : lecture 
de disque dur, contrôle de positionnement dans l’espace, codage de position 
et d'angle, mesure de courant électrique sans contact, imagerie magnétique, 
contrôle non destructif, imagerie médicale, ils permettent en effet de détecter 
et mesurer un champ magnétique créé par un objet magnétique (aimant ou 
bobine) ou par un courant électrique (à travers son champ magnétique 


rayonné). 


Il y à environ une quinzaine de types de capteurs magnétiques. Certains sont 
des capteurs de champ (c’est-à-dire des capteurs mesurant directement la 
valeur du champ magnétique selon un ou plusieurs axes), d’autres, des 
capteurs de flux (c’est-à-dire mesurant l’intégrale du champ passant à travers 


une surface). 


Actuellement, la plupart des entreprises dominant le marché se sont lancées 
dans ce développement, et il est probable que dans un proche avenir les 
capteurs magnétiques basés sur l’électronique de spin auront une place 


prépondérante dans cet ensemble. 
I.5.3 Mémoire magnétique à accès aléatoire (MRAM) 


De nombreuses recherches ont été menées pour trouver une mémoire à l’état 
solide qui pourrait concurrencer les mémoires mécaniques telles que les 


disques durs. 


Les mémoires magnétiques à accès aléatoire sont un nouveau type de 
mémoire informatique à priori très rapides non volatiles et avec une 
endurance quasi infinie. Dans cette mémoire, les données sont stockées sous 
forme d’aimantation sur des nœuds d’un réseau plans bidimensionnel à Pétat 
solide (une architecture matricielle). Comme nous allons le voir, ces lignes 
servent à la fois à lire le bit et à l'écrire. La figure I.6, représente une vue 


schématique de cette matrice mémoire. 
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MT] element 


Figure I.6 : Architecture cross point pour une mémoire MRAM 


constituée d’une matrice de jonctions tunnels magnétiques [16]. 


Le principe de fonctionnement du MRAM est simple, comme on le voit dans 
la figure L.7, il y a un maillage de lignes électriques où à chaque intersection se 
trouve une jonction magnétique tunnel qui forme un bit. Ce bit est lu par effet 
TMR et il est écrit par le champ magnétique induit par un courant placé à 
proximité sur le "bit line". Plus la TMR de la jonction tunnel magnétique est 
importante, plus les deux niveaux de résistance seront séparés, et donc mieux 


les bits "0" ou "1" seront définis. 


transistor 


Figure [1.7 : MRAM 4 Mb commercialisée par Freescale depuis juillet 2006 


avec vue agrandie d’une cellule simple [17]. 


En pratique, une TMR de l'ordre de 100 % assure la bonne discrimination de 
ces niveaux et ces mémoires, en cours de développement, combine de 
nombreux avantages potentiels : vitesse, endurance et non volatilité. Ce 
concept de mémoire, mis sur le marché en 2006 par Freescale, est en plein 


Essof. 
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L.6 Les semi-conducteurs magnétiques dilués DMS comme une 


nouvelle classe des matériaux pour l’électronique de Spin 


Après la découverte de la GMR et de la TMR [6,7], la conception et la 
fabrication de nouveaux matériaux ont été accélérées et ont fait l’objet de 
beaucoup de recherches aussi bien sur le plan théorique (fondamental) que 


sur le plan application. 


Néanmoins, la problématique de linjection d’électrons avec une direction de 
spin contrôlée dans l’électronique traditionnelle à poussé au développement 
futur de la discipline à travers la fabrication de structures et de composants 
hybrides, associant métaux magnétiques et semi-conducteurs. En effet, depuis 
près de quatre décennies, les méthodes connues pour générer un déséquilibre 
entre les populations de porteurs de spin minoritaire et majoritaire dans un 


semi-conducteur sont l'injection par pompage optique et l'injection électrique. 


La technique de pompage optique utilise une lumière polarisée circulairement 
pour exciter les électrons depuis la bande de valence vers la bande de 
conduction, c’est une méthode de polarisation des porteurs dans un semi- 
conducteur, antérieure aux développements récents de la spintronique. 
Cependant, l'injection optique a rencontré un obstacle fondamental pour 


créer des dispositifs efficaces industrialisables. 


La technique de l'injection électrique se présente comme la solution la plus 
appropriée dans laquelle on utilise le courant d’un matériau fortement polarisé 
dans un semi-conducteur à l’aide d’une jonction métal ferromagnétique/semi- 
conducteur (MF/SC), et le résultat souhaité est que ce courant doit conserver 
sa polarisation lorsqu'il traverse le semi-conducteur; c’est le cas des composés 
légèrement dopés par des ions d’éléments de transition. Dans ce contexte, les 
semi-conducteurs magnétiques dilués (DMS) [17,18] répondent bien à ce 


besoin et possèdent les caractéristiques souhaitées. 
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Les recherches actuelles sur les semi-conducteurs magnétiques dilués portent 


sut deux axes : 


" D'une part, une grande valeur de la température de Curie pour 
améliorer la qualité d’un matériau donné, ou en cherchant carrément 
des matériaux plus favorables pour permettre leur utilisation dans des 


conditions normales de température (température ambiante). 


" D'autre part, fabriquer des composants élémentaires de l’électronique 


de spin à partir des semi-conducteurs magnétiques dilués actuels. 


En effet, le comportement semi-conducteur dans les matériaux 
ferromagnétiques a été découvert dans les années 60. On peut citer les 
chalcogénures d’europium [19-21] comme le EuO (Te = 69 K) et Eus (Te — 
16,5 K) et les semi-conducteurs de type spinelle [22-25] comme le CdCrS4 
(Tc=84 K), HgCnSes (Tc=106K) et CdCrSes (Tc=130 K). 


Ces matériaux sont caractérisés par une forte concentration d'atomes 
magnétiques et sont appelés CMS (Concentrated Magnetic Semiconductors) ; 
ces éléments magnétiques forment un réseau périodique figure L.8 (A), donc 
une structure cristalline très différente de celle des semi-conducteurs « 
traditionnels » tels que Si et GaAs ; de plus, leur température de Curie est très 
basse, inférieure à 200 K , et par conséquent, leurs applications dans les 


technologies semi-conductrices modernes sont limitées. 


Pour pallier à cette contrainte, dans les années 70 et au début des années 80, 
le groupe de Robert Galazka [26,27] à réalisé les premiers DMS, plus proches 
que possible des semi-conducteurs usuels de la microélectronique par 
l'introduction du manganèse (Mn), comme impureté magnétique, dans la 
matrice d’un semi-conducteur hôte de la famille II-VI. Dans ce type de 
matériaux, l'élément magnétique vient en substitution du cation de l’hôte 
semi-conducteur figure [.8 (B), et présente des propriétés magnétiques très 


couplées aux propriétés électroniques. 
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Le premier composé étudié dans la famille des semi-conducteurs de type III- 
V fut l’arséniure d’indium (InAs), dopé au Mn sous forme de couche mince 
[28]. Le groupe de H.Ohno à rapporté en 1989 l’existence d’une phase 
homogène d’Ini,Mn,As ferromagnétique [28], puis ils ont montré en 1992 que 


le ferromagnétisme était induit par les trous [29]. 


Les travaux de recherche publiés sur ces matériaux ont conduit à d’autres 
recherches sur les semi-conducteurs du groupe III-V mais la température de 


curie reste toujours la limitation majeure de ces DMS. 


Figure I.8 : Classifications des semi-conducteurs : (A) Semi-conducteurs 
magnétiques concentrés (CMS) : les éléments magnétiques forment un sous 
réseau de semi-conducteurs. (B) Semi-conducteurs magnétiques dilués : les 
atomes magnétiques sont aléatoirement positionnés dans le réseau du semi- 
conducteur hôte. (C) Semi-conducteurs non magnétiques, dans ce cas, le 


matériau est un semi-conducteur binaire (d’après Ohno) [28,29]. 


Nous pouvons distinguer différents types de semi-conducteurs magnétiques 
dilués selon la matrice semi-conductrice hôte avant le dopage par une 
impureté magnétique ; et comme nous l'avons dit précédemment, les DMS à 
base des semi-conducteurs II-VI [27] étaient les premiers étudiés par rapport 


aux DMS des groupes II-V, IV et II-VI-VI2 récemment étudiés. 


La famille des semi-conducteurs magnétiques dilués de la colonne IV est 
généralement compatible avec les technologies actuelles de la 


microélectronique mais jusqu’à présent elle à fait objet de peu d’études. Ces 
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matériaux sont composés d’un seul type d’atomes (Si et Ge) dopés par 
exemple au Cr, Mn, Ni. On peut noter que les premiers résultats concluants 
sur des films de GeMn remontent à 2002 [30] ; ces résultats ont montré une 
relation directe entre la concentration en Mn et la température de Curie (qui 
varie entre 25 K et 116 K) dans des films minces déposés par épitaxie par jet 


moléculaire. 


Les semi-conducteurs magnétiques dilués de la famille IV-VI sont des semi- 
conducteurs composites fabriqués à partir d'un ou de plusieurs éléments de la 
colonne IV du tableau périodique des éléments et d'un ou de plusieurs 
éléments de la colonne VI, ils présentent une très grande densité de porteurs 
suffisante pour que les interactions ferromagnétiques qu’ils véhiculent soient 
plus fortes que l’antiferromagnétisme. Toutefois, ces matériaux ne sont pas 
compatibles avec ceux de lélectronique moderne, car leur structure 


cristallographique est trop différente. 


La famille des DMS à base des semi-conducteurs II-VI-VI (appelés aussi les 
chalcopyrites) à l'avantage de donner des températures de Curie supérieures à 
300 K, mais ces composés sont difficiles à incorporer dans les structures IV 
ou III-V et ne sont pas de bons candidats pour la réalisation de dispositifs en 
raison de la difficulté de leur élaboration. Cependant, seulement des résultats 
pour les matériaux massifs (ZniMn,;GeP: et Cdi:Mn;GeP2) ont été 


rapportés (Ie 310 K) [31,32]. 


Comme on l’a vu, les premiers travaux ont porté sur les semi-conducteurs II- 


VI, la recherche s’est ensuite orientée vers les semi-conducteuts IIT-V. 


En effet, dans les semi-conducteurs II-VI le dopage avec les métaux de 
transition a montré de bonnes propriétés de solubilité, à cause de la valence 
identique entre les éléments de la colonne II et les métaux de transition où, 
en effet, la solubilité atteint les 80 % [33]. Ces mêmes valeurs de solubilité ont 


même été prédites pour le manganèse dans les composés ZnT'e et CdTe (Mn 
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remplaçant alors Zn et Cd), comparativement aux semi-conducteurs III-V 
(GaAs, InAs, GaN, AIN....), dans lesquels les valeurs de solubilité des métaux 
de transition sont beaucoup plus faibles, ce qui donne des difficultés à obtenir 
des matériaux suffisamment concentrés pour avoir des propriétés 
magnétiques intéressantes. Cependant, même si ces prédictions théoriques 
ont été confirmées expérimentalement, les propriétés magnétiques de ces 


matériaux étaient en deca des aspirations des chercheurs. 


Néanmoins, ils possèdent de nombreux avantages : ils possèdent des 
structures cristallographiques compatibles avec des dispositifs existants 
comme le GaAs, et ils sont efficaces dans la polarisation et l'injection de spins, 
ils sont caractérisés par des durées de vie plus importantes (jusqu’à 100 ns) 
des porteurs polarisés permettant leur transport sur plusieurs centaines de 
nanomètres. Bien évidemment, ils ont des propriétés magnétiques et 


présentent une phase ferromagnétique au-delà de la température ambiante. 
L.6.1 Les applications des semi-conducteurs magnétiques dilués 


L'utilisation des semi-conducteurs dans l’électronique de spin a favorisé une 
intense recherche qui donnera certainement une nouvelle génération de 
composants électroniques. Tout d’abord, un des avantages de la spintronique 
(dans les semi-conducteurs) est le temps de vie du spin qui est plus long que 
dans les métaux, et donc l'information de spin peut se propager sur des 


distances plus grandes. 


Ainsi, l'intégration du spin dans les semi-conducteurs favorisera un vaste 
champ de dispositifs très efficaces, et constitue un nouveau défi. Parmi ces 
dispositifs les plus célèbres le Transistor à effet de Champ de Spin (Spin-FET) 
et le transistor de Spin à base de Diodes électroluminescentes « Spin Light 


Emitting Diode » (Spin-LED). 
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L.6.1.1 Le transistor à effet de champ de spin (Spin-FET) 


Le transistor à effet de champ (Spin-FET) est un futur dispositif de 
spintronique semi-conducteur doté de performances meilleures par rapport à 


celle obtenue avec la technologie actuelle des transistors. 


Un schéma illustratif du Spin-FET est représenté sur la figure 1.9. Ce 
transistor est composé de deux contacts ferromagnétiques (source et drain), 
qui prennent en sandwich la région semi-conductrice, les contacts 
ferromagnétiques contiennent principalement des électrons polarisés en spin 
et jouent le rôle de polariseur et d'analyseur comme décrit en 1990 par les 


deux chercheurs américains Datta et Das [34]. 


Le champ électrique induit 
une rotation de spins 


Courant bloqué 
à l'interface 


Figure L.9 : Schéma de fonctionnement d’un transistor Spin-FET. 


Le contact de source ferromagnétique injecte des électrons polarisés en spin 
dans la région semi-conductrice (Si, GaAs, etc). En raison de l'interaction 
spin-orbite non nulle, le spin de l'électron se propage dans le canal. Au contact 
du drain, seulement les électrons avec le spin aligné à l'aimantation du drain 
peuvent facilement quitter le canal et contribuer au courant. Ainsi, le courant 
total à travers le dispositif dépend de l'angle relatif entre la direction 


d'aimantation du contact de drain et la polarisation de spin électronique à 
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l'extrémité du canal semi-conducteur. Une tension de grille crée dans le 
canal un champ électrique qui va provoquer ou non la rotation des 
spins électroniques par le phénomène de précession de Rashba. La 
modulation du courant est obtenue en accordant la force de l'interaction spin- 


otbite dans la région semi-conductrice et donc le degré de précession de spin. 


I.6.1.2 Le transistor de Spin à base de Diode électroluminescente « Spin 


Light Emitting Diode » (Spin-LED) 


Les sources lumineuses à polarisation de spin constituent une nouvelle classe 
de dispositifs dans lesquels la recombinaison radiative de porteurs polarisés 
en spin entraîne une luminescence présentant une polarisation circulaire nette. 
Le transistor de Spin à base de Diode électroluminescente est un composant 
électronique qui permet de produire de la lumière quand il est parcouru par 
un courant électrique, il est constitué d’un semi-conducteur photoémissif 
dont les caractéristiques électriques correspondent à celles d’une diode. La 
mesure de la polarisation circulaire de la lumière émise par la spin-LED 
renseigne directement sur la polarisation en spin des électrons injectés dans 


le semi-conducteur. 


Dans un tel dispositif (figure 1.10), la première couche de GaAs (substrat) est 
dopée n, sur laquelle est déposée par épitaxie une autre couche de GaAs dopée 
n, les trois couches suivantes GaAs-InGaAs-GaAs ne sont pas dopées, la 


dernière couche est constituée du DMS GaAs : Mn dopé p. 
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E] 


NO > ® © 


EL intensity (arbitrary units) 
Polarization (%) 


125 1.30 1.35 1.40 1.45 1.50 1.55 1.60 
Energy (eV) 


Figure I.10 : Fonctionnement d’une diode Spin-LED basée sur le GaAs. (a) 
Schéma d’une Spin-LED basée sur du GaAs multicouche. (b) Courbe de 
l'intensité électro lumineuse émise par la Spin-LED avec T = 6 K et H = 1000 


Oe, la distance d sur le schéma vaut 20 nm [35]. 


Notons que le courant de spin-polarisé est transporté dans la couche non 
magnétique GaÂs au puit quantique InGaÂs où les trous polarisés se 
recombinent avec les électrons de même spin provoquant l'émission de 


photons polarisés circulairement. 


La production en grande série de Spin-LED vient seulement de commencer, 
par exemple dans la référence [36], l’injecteur CoFe/MeO était déposé sur 
une diode photoluminescente (LED) à base de puits quantiques 
GaAs/AIGaAs/GaAs et aussi d’autres types de dispositifs ont aussi été 
réalisés avec un semi-conducteur II-VI tel que ZnMnSe sous champ 


magnétique intense et basse température pour aligner les spins des ions Mn 
(371: 


[1.7 Les demi-métaux dans les DMS et leur intérêt pour l'électronique 


de spin. 


Récemment, plusieurs nouveaux matériaux magnétiques ont été découverts 


comme les demi-métaux ferromagnétiques. Ce caractère a été présenté pour 
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la première fois par de Groot et al [38] en 1983 quand ils ont étudié la structure 


de bande des alliages semi-Heusler NiMnSb et PtMnSb. 


La propriété exceptionnelle d’un demi-métal est qu’il est semi-conducteur 
pour une orientation du spin et métal pour la seconde. Il a une polarisation 


de spin égale à 100 % au niveau de fermi (Er) (P=1), voir figure 1.11. 


A. Non ferromagnetique B. Ferromagnétique C. Demi-métallique 


N(E) N(E) 


P-=0% P<50% 

Fe, Co, Ni 

Figure I.11 : (A) Représentation schématique des densités d'états et 
£ P q 

polarisation en spin d’un métal non ferromagnétique, (B) d’un matériau 


ferromagnétique et (C) d’un matériau-demi métallique. 
La polarisation en spin d'un système est définie par l’équation (1.5) 


_ N'(Er) — N°(Er) (LS) 


PME) + NE) 


N'(E-) et N!(Ek) sont les densités d’états au niveau de fermi pour les deux 


directions de spins up et down respectivement. 


Les demi-métaux ne doivent pas être confondus avec les ferromagnétiques 
forts comme le Co ou le Ni. En effet, les bandes ‘3d° du Co ou Ni sont bien 
polarisées en spin à 100 % mais les bandes ‘4s’, qui se trouvent au niveau de 
Fermi, ne sont pas polarisées ; des électrons « up » ou « down » sont donc 
présents au niveau de Fermi. Pour obtenir un demi-métal, il faut alors une 


hybridation des bandes ‘34° et ‘4s’ pour que le niveau de Fermi ne se trouve 


27 


CHAPITRE I. LE ROLE DU SPIN DANS L’ELECTRONIQUE 


plus dans la bande ‘4s’. C’est pour cela qu'aucun matériau constitué d’un seul 


atome n’est demi-métallique. 


Depuis l'étude de Groot et al [38], d’autres demi-métaux ont été prédits dans 
d'autres structures, tels que le graphène [39], les oxydes métalliques Fe3O4 [40], 
CrO2 [41, 42], les demi-Heusler Co2MnSi [43], Co2FeSi [44], les alliages 
pérovskites LaorStoaMnO; [45], SrFeMoO« [46] et différents semi- 
conducteurs binaires dopés avec des impuretés magnétiques dans la phase 


rock-salt [47] et aussi dans la phase zinc-blende [48]. 


D'une manière plus claire, les demi-métaux ferromagnétiques peuvent 
répondre aux trois grands défis technologiques dans le domaine de 


l'électronique de spin. 


On commence avec le défi de la magnétorésistance (MR) qui est un effet 
principal dans la spintronique et concerne la variation de la résistance 
électrique d’un conducteur sous Papplication d’un champ magnétique. 
Aujourd’hui, des valeurs très élevées de MR sont hautement recherchées pour 
le développement des dispositifs spintroniques tels que les mémoires 
magnétiques (MRAM). Il existe plusieurs types de MR mais les plus utilisées 
sont la magnétorésistance géante (GMR) et la magnétorésistance tunnel 
(TMR). L'effet GMR, observé dans les systèmes de multicouches qui sont 
composées alternativement d’une couche magnétique et d’une couche 
métallique non-magnétique, est utilisé généralement dans les capteurs 
magnétiques et dans les têtes de lecture des disques durs. Un autre effet de 
magnétorésistance similaire à la GMR, est la TMR qui a été observé dans les 
structures type jonctions magnétiques tunnel (MT] : magnetic tunnel 
junctions) obtenues en remplaçant la couche séparatrice dans les structures 
GMR par une fine couche isolante (barrière) comme ALO3 , MeO. Ce taux 
de MR est proportionnel à la polarisation en spin, P, des électrodes selon le 


modèle de Jullière d’où la nécessite d’utiliser des électrodes demi métalliques. 
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L'application principale de la TMR se situe dans la mémoire magnétique à 


accès aléatoire (MRAM : magnetic random access memory). 


L'autre défi de ce domaine concerne l’injection d’un courant polarisé en spin 
d’un matériau ferromagnétique vers un autre matériau semi-conducteur où la 
désadaptation de résistance entre ces deux matériaux constitue l’obstacle 
majeur pour cette injection. Il existe deux solutions pour obtenir une injection 
efficace de spin dans les semi-conducteurs. La première est d’utiliser un semi- 
conducteur magnétique à la place d’un ferromagnétique. Cependant, la 
température de ces semi-conducteurs ferromagnétiques comme GaMnAs 
reste largement au-dessous de la température du Curie. La deuxième solution 
consiste à utiliser un matériau ferromagnétique ayant une forte polarisation 


en spin à la température ambiante tel qu’un demi-métal. 


Le troisième et dernier défi important est le transfert de spin (STT : spin 
transfer torque) permettant le contrôle et la manipulation de l’aimantation 
d’un matériau ferromagnétique nanométrique par un courant polarisé en spin. 
Depuis la découverte du phénomène STT en 1996, un intérêt considérable a 
été prêté au développement des dispositifs à base de STT tels que les MRAM 
ou les oscillateurs hyperfréquences. En effet, écrire une information dans les 
mémoires magnétiques en utilisant un courant polarisé en spin simplifiera leur 
architecture et réduira leur consommation électrique. Cependant, pour être 
plus compétitif et pour augmenter plus significativement la densité de 
stockage des mémoires et la vitesse de retournement de laimantation, la 
densité de courant doit être réduite d’avantage. Cette densité de courant est 
proportionnelle au coefficient d'amortissement de Gilbert et inversement 
proportionnelle à la polarisation en spin des matériaux utilisés. L'utilisation de 
ceftains  semi-conducteurs magnétiques  dilués pourrait réduire 


significativement cette densité de courant. 
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L.8 Les interactions entre les porteurs localisés et délocalisés dans les 


DMS 


Pour déterminer l’origine du ferromagnétisme dans les semi-conducteurs 
magnétiques dilués, la compréhension des interactions magnétiques est bien 
évidement nécessaire. Ces interactions sont principalement dues à la 
coexistence de deux aspects électroniques : une contenant des électrons 
présents dans la bande de valence ou dans la bande de conduction, et des 


électrons d'impuretés magnétiques portant un moment magnétique. 


Parmi les interactions qui dominent les propriétés magnétiques des semi- 
conducteurs magnétiques dilués, nous notons : les interactions entre les 


électrons de type ‘d’ des ions magnétiques et les porteurs délocalisés. 


Les interactions ‘s-d’ et ‘p-d’ viennent de la présence de deux systèmes 
électroniques, un système de spin lié aux électrons ‘s-p’ des bandes de valence 
et de conduction du semi-conducteur, l’autre système lié aux électrons ‘3d? 


des ions magnétiques de l’élément de transition. 


L’Hamiltonien de Heisenberg qui décrit l’échange ‘sp-d” est le suivant [49] : 


Hsp-a 7 dis (r —R;)S;o (L.6) 


" S;est le spin des ions magnétiques (en position R;). 
" œest le spin des porteurs délocalisés (en position r). 


" Jp-a: l'interaction d’échange entre les porteurs délocalisés et les 


électrons des ions magnétiques. 


Dans l'hypothèse où les ions magnétiques sont répartis aléatoirement dans la 
matrice semi-conductrice, nous pouvons remplacer le spin S; par sa 


moyenne (S), c’est l’approximation du champ moyen. Donc, l’hamiltonien 
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d'échange entre les ions magnétiques et les porteurs délocalisés dans la bande 


de conduction et dans la bande de valence peut s’écrire comme suit : 


Ha = —XNoa(S)o (L.7) 
Ha = —XNoB(S)o (L.8) 


Où : 
" N est la concentration volumique des sites cationiques. 
" a l’intégrale de recouvrement entre les orbitales ‘s’ et ‘d’. 
" B l'intégrale de recouvrement des orbitales ‘p’ et ‘d’. 


[9 L'état de l’art sur les semi-conducteurs magnétiques dilués ILA-VIs 


à base de MgSe. 


La recherche des matériaux multifonctions, à propriétés étendues pour une 
application variée avec une efficacité optimale, fût l'objectif de la 
communauté scientifique. Pour améliorer l'efficacité des différents dispositifs 
magnéto résistifs, de nombreuses recherches ont été menées ces dernières 
années sur les matériaux ferromagnétiques, et particulièrement sur les semi- 


conducteurs magnétiques dilués. 


Depuis des années, les composés Ils — Vs et IA — VIs dopés par un élément 
de transition ont été largement étudiés. Les propriétés magnétiques très 
particulières, dues à l'addition d’impuretés, apparaissent du fait que les états 
‘3d d’impuretés sont partiellement remplis. Les états ‘s-p” des binaires 
IB — Vis et Ia — Vis s'hybrident fortement avec les états ‘3d’ de l'ion 
magnétique pour donner naissance à des propriétés électroniques et 


magnétiques intéressantes au plan technologique. 


Dans ce contexte, les alcalo-chalcogénures des groupes IA-VIs (IA = Be, Mo, 
Ca, Sr et Ba, VI = O,S, Se et Te), ont suscité un intérêt technologique 


considérable, en raison de leurs propriétés physiques importantes ; ils sont 
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caractérisés par une large bande interdite (Gap) et présentent également les 
propriétés des semi-conducteurs, ils sont aussi sélectionnés comme matériaux 
prometteurs dans des dispositifs électroniques et  optoélectroniques 


fonctionnant sur de larges régions spectrales [50]. 


En particulier, les chalcogénures de magnésium MeX (X5S, Se, et Te), sont 
un bon choix pour les revêtements protecteurs en raison de leur dureté, de 
leur point de fusion élevé, de leur conductivité thermique élevée et de leur 


module de compressibilité [51-53]. 


Plus précisément, le sélénium de magnésium (MgSe) est un semi-conducteur 
binaire intéressant de la famille ITA-VIs, ses applications vont de la catalyse à 
la microélectronique [51]. Il peut être utilisé à haute température, dans les 
diodes laser (LD) et les diodes électroluminescentes (DEL) à base de ZnSe en 
tant que matériau de revêtement [51,52], et comme dopant de type n 


et p [54-561. 


Bien que, le composé MgSe présente plusieurs phases cristallines [53] : 
rocksalt NaCI (B1), zinc blende (B3), wurtzite (B4), et nickel arsenide NiAs 
(BB), il a fait au cours des dernières années, l’objet de nombreuses études 


expérimentales et théoriques. 


Sur la base des caractéristiques des composés susmentionnés, nous avons 
essayé de collecter les études disponibles dans la littérature qui a porté sur les 
effets du dopage par des métaux de transition (TM) sur les semi-conducteurs 
IlAa-Vis à base de MoSe, notant qu'il n'y à pas encore de résultats 


expérimentaux suf ce COMpPOosé. 


MéSajjad et al [57], à étudié théoriquement avec le code WIEN2K les 
propriétés structurales des composés Mg: V,Se/Te (x= 0.25, 0.50, 0.75) dans 
la phase zinc blende en utilisant l’approximation du gradient généralisé 
(GGA) de Wu-Cohen. Ce groupe de chercheurs ont confirmé le 


ferromagnétisme demi-métallique (HMP) dans le MgSe dopé au Vanadium 
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avec différentes concentrations en utilisant le potentiel modifié de Becke et 


Johnson (TB-mBjJ). 


Noor et al [58], via l’approche de Wu-Cohen GGA, ont exploré le caractère 
HMF dans le composé MseSe dopé au Cr et ont montré par l'optimisation de 
la variation de l’énergie en fonction du volume que ces alliages sont stables 
dans la phase ferromagnétique. D’autres calculs de « premier principe » 
réalisés avec l’approximation du gradient généralisé (GGA) de Perdew-Burke- 


Ernzerhof (PBE) ont mené à des résultats similaires. 


Cependant, dans le cas de MoSe dopé au Ti pour une concentration fixe 
(x=0.25), la stabilité magnétique de ce composé à été vérifiée par la différence 
d'énergie totale (AE) de l’état ferromagnétique à l’énergie totale de l’état verre 
de spin, où la valeur négative de (AE) indique que la phase FM est stable. En 
outre, les propriétés élastiques réalisés par Q. Mahmood et al [59] ont montré 
que ce composé a un fort comportement ductile, en plus de l'existence du 


comportement demi-métallique. 


W.Tanverr et al [60], ont étudié l’effet du dopage par différents éléments 
magnétiques tels que Fe, Co, et Ni avec une concentration fixe (x=0.25) sur 
les propriétés physiques des composés MgS/Se/Te. Les résultats des 
propriétés électroniques ont conclu que le MgSe dopé au Fe et au Ni présente 
un comportement  demi-métallique tandis qu’un comportement 
ferromagnétique semi-conducteur a été obtenu pour le MgSe dopé au Co. Le 
spectre optique a montré que la lumière maximale absorbée se situe dans la 
région visible et dans la région UV du spectre électromagnétique. Cette 
caractéristique optique fait de ces composés MeXSe (X= Fe, Co, et Ni) des 
candidats potentiels pour les applications dans des dispositifs 


optoélectroniques et spintronique. 
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CHAPITRE II. CONTEXTE THEORIQUE DE LA 
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II.1 Mécanique quantique de la physique de la matière condensée 


La détermination de la structure électronique d'un matériau passe 
nécessairement par la résolution de l’équation de Schrôdinger (IL.1) [1] d’un 
système à N corps dans lequel de multiples interactions peuvent avoir lieu : 
les interactions entre électrons, les interactions entre les électrons et les 
noyaux et les interactions entre noyaux. Prendre en considération toutes ces 
interactions rend impossible la résolution de l’équation de Schrôdinger à 
l'exception du cas particulier de systèmes ne contenant qu’un seul électron et 


un seul noyau tel que atome d’hydrogène ou la molécule H7. 
HD = Ep (1) 
Où 


= est l'état propre (fonction d'onde totale de l'état fondamental à 


déterminer). 
= Hest l'hamiltonien total du système. 
“ Eest la valeur propre de l'énergie de l'état propre. 


L'expression de l’'Hamiltonien H (en unités atomiques) est donnée par 


équation (IL.2). 
2 1 V2 1 1 

n- on ns 

ZM 24m GR 24-56 


y aa 
24 \R;, —R; (IL2) 


Où 


" Met m sont respectivement les masses du noyau et de l’électron. 
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" Zest la charge nucléaire. 


" Ret f sont respectivement les vecteurs de position du noyau et de 


l'électron. 


Les cinq termes de l’équation (IL.2) s’écrivent : 


nr IR À; (1.3) 


; PR RS (11.4 
14 =-) ——:0 = Ù — 
dl RE 24 ñ-f] 


Donc l'hamiltonien peut être écrit sous la forme suivante : 
À =T+T, + Von + ne + Üee (LS) 


Les premier et deuxième termes de (IL.5) sont les énergies cinétiques de 
l'électron et des noyaux, respectivement, les troisième et quatrième termes 
sont respectivement la répulsion noyaux-noyaux et l’attraction coulombienne 


noyaux-électrons, et le cinquième terme est la répulsion électron-électron. 


Comme mentionné auparavant, l'équation de Schrôdinger est trop difficile à 
résoudre avec cet hamiltonien et des approximations s'imposent pour 


résoudre les systèmes complexes. 
II.2 L’approximation de Born-Oppenheimer 


La première approximation est celle de Born Oppenheimer [2] qui tient 
compte des différences des vitesses de déplacements des électrons et des 
noyaux pour raison de leurs masses. Tenant compte du fait que la masse des 
électrons est relativement trop faible devant celle des noyaux, il est par 
conséquent possible de négliger l'énergie cinétique des noyaux devant celle 


des électrons. 
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Dans ces conditions T, est pris égal à zéro et V,, est une constante. Le terme 
Pre s'écrit alors V4 car il s'agit maintenant d'un potentiel «externe» dû aux 


noyaux stationnaires auxquels sont soumis les électrons. 


Suite à ces nouvelles considérations, l'Hamiltonien complet prend la forme 


suivante : 
B =T + Vo + Ue (IL.6) 


Cet Hamiltonien est mieux allégé que l'original sauf que le problème du grand 
nombre de cofps reste toujours posé d’où la persistance de la difficulté à 
résoudre l’équation (IL.6). D'autres approximations non traitées ici, telle que 
celle de Hartree-Fock [3,4], ont tenté de trouver une solution à ce problème 


mais vainement. 
IL.3 Approximation de Hartree et de Hartree-Fock 


Appelée aussi Papproximation mono électronique [3,4], elle consiste à 
considérer chaque électron comme étant indépendant mais subissant un 
potentiel effectif, Ve qui doit tenir compte de toutes les interactions 
électroniques. Dans l’approximation de Hartree, la fonction d’onde du 
système à n électrons se ramène à un produit de n fonctions d'ondes à un 


seul électron, qui sont solutions de Péquation suivante : 


1 
|-574 + Vers (M) D(r) — E;@(r) (LT 


Le potentiel V,,-(n) est donné par la somme de l'interaction électron-noyaux, 
Pe et d’un potentiel Üz, dit de Hartree (ou champs moyen), qui doit tenir 
compte des interactions coulombiennes répulsives entre un électron de 
cordonnée 7% plongé dans le champ moyen des autres électrons de 


cordonnées 7 . 


1 ” re 
[578 + (À) + Ouf T)]|0 = 500 9 


41 


CHAPITRE IT. CONTEXTE THEORIQUE DE LA MATIERE CONDENSEE 


En 1930, Fock a noté que la fonction d’onde de Hartree viole le principe 
d'exclusion de Pauli, car il néglige un terme important qui est le terme 
d'échange ce qui suppose que les électrons sont considérés sans spin, et les 


solutions ne sont pas antisymétriques. 


Fock a remplacé alors la fonction d’onde par un déterminant de Slater de 
fonctions d’onde mono électroniques, qui sont antisymétriques dans 
l'échange de deux électrons et le principe de Pauli ici est donc satisfait. Les 
équations de Hartree-Fock contiennent, en plus des termes de Hartree, un 


terme dit d'échange, c’est lui qui fait la différence entre les deux méthodes. 


1 es LS 
(PE + Mel À) + Ou) + xD] EG = 600 ae 


L’approximation de HF est assez précise pour décrire les atomes et les 
molécules mais moins appropriée pour le système qui comprend un grand 
nombre d'électrons ce qui est le cas dans les solides. Une théorie plus efficace 


pour les solides est la théorie de la fonctionnelle de la densité. 
II.4 La théorie de la fonctionnelle de la densité (DFT) 


La DFT est une théorie mise en œuvre pour apporter des solutions au 
problème à plusieurs corps tout en contournant la solution qui passe par 
l'équation de Schrüdinger basée sur un Hamiltonien dans sa forme la plus 
complète. Cette théorie a comme fondements les deux théorèmes (Postulats) 


de Hohenberg et de Kohn, établis en 1964 [5]. 


Théorème I : I] existe une relation un-à-un unique entre le potentiel externe Vore (À) ét 


la densité électronique p(F). 


Enoncé : La valeur d'espérance fondamentale de tout élément observable, 
telle que À est une fonction unique de la densité d'électrons à l'état 
fondamental. Cela indique que la densité de l'état fondamental contient autant 


d'informations que les fonctions d'onde de l'état fondamental. 
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Théorème II : L'énervie totale de l'état fondamental E\p] est une fonction de la densité 


d'état fondamental p(r) et est de la forme : 
Elp] = (ÏT + OelD)HD|Pex|D) (L.10) 


Elo] = F[p] + Î PCF) Vexe(F)d? t-11) 


Où F[p] est une fonction universelle pour tout système à plusieurs électrons, 
Ep] est variationnelle au sens qu’elle est minimisée par la densité de Pétat 


fondamental correspondant à V,. 


Enoncé : la fonctionnelle F[p] est indépendante de toute information sur les 


noyaux et est donc universelle pour tout système à plusieurs électrons. 


En résumé, ces théorèmes stipulent que le potentiel de tout système est 
uniquement défini par la densité de l’état fondamental, et l'énergie minimale 
totale de l'état fondamental correspond à la densité de l'état fondamental. 
Nous avons seulement besoin de résoudre la densité car elle contient autant 
d'informations que la fonction d'onde. Nous avons aussi besoin que de la 
fonctionnalité universelle F[o] pour résoudre tout le problème ; 
malheureusement, une expression explicite pour F[p] n'est pas connue. Kohn 
et Sham ont fourni une méthode pour trouver la densité de l’état fondamental 


pour pouvoir résoudre l'équation (IL.9). 
IL.5 Équations de Kohn-Sham 


L'idée de Kohn et Sham [6] était de cartographier le véritable système d'auto- 
interaction pour un système auxiliaire de particules sans interaction. Cela 
permet d'inclure toutes les interactions dans un seul terme et simplifie les 


calculs. Les fonctions universelles peuvent être séparées en trois termes 


F[p] = To + Vy + Ve (I.12) 
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"  Tiest l'énergie cinétique du gaz d'électrons sans interaction. 
" V; est le potentiel de Hartree qui est dû à l'interaction classique de 


Coulomb entre les électrons. 


" est le potentiel d'échange-corrélation défini comme suite : 


6Exe (1.13) 


Les deux premiers termes de (11.12) sont connus exactement. Le terme 
échange-corrélation contient la différence entre les énergies cinétiques exactes 
et non interactives ainsi que les effets non classiques d'échange et de 
corrélation. C'est ce terme qui inclut les interactions complexes entre 
électrons et peut prendre différentes formes qui vont définir les 


approximations que nous allons voir plus loin. 
L'hamiltonien d’un système s’écrit alors : 
Bkslp] = To + Va + ec + Voxt (L.14) 


Le système auxiliaire est constitué de particules individuelles n'interagissant 
pas, de sorte que la fonction d'onde globale de l'état fondamental de ce 
système peut être donnée avec exactitude par un ensemble de fonctions 


d'onde à une seule particule. 


La densité d’un tel système est construite à partir des fonctions d'onde ;(r) 


d'une particule comme suit : 


A (II.15) 
pE) = > lei PIE 


Nous pouvons, par conséquent, résoudre les équations de Schrôdinger à une 


seule particule d’une manière simple une fois que H4s{p] est connu: 


AxsDi — Ep; (1.16) 


44 


CHAPITRE IT. CONTEXTE THEORIQUE DE LA MATIERE CONDENSEE 


Les équations (11.14), (IL.15), et (11.16) sont les équations de Kohn-Sham qui 
peuvent être résolues de manière itérative. D'abord, une première estimation 
pour la densitépyest faite et mise dans l'équation  (IL.14). 
L’Hamiltonien Hxsest formé et mis en équation (IL16) pour donner 
l'ensemble des g;. Finalement, les y; sont utilisées dans (11.15) pour donner 


une nouvelle densité p:. 


Le cycle se poursuit jusqu'à ce qu'une densité auto cohérente finale soit 
trouvée. Cette densité finale est dite cohérente avec l'Hamiltonien et 
cofrespond à la densité de l'état fondamental du système réel. Un schéma 


décrivant les principes de la DFT est donné par la figure IL.1. 


Calcul du potentiel effectif 
Va = Veste + Virfn] + Vxc[n] 


Résolution des équations de Kohn-Sham 


[9 +16 ét = tnt) 


Calcul de la nouvelle densité électronique 


n°(r)= XF 


Figure IL.1 : Les principes de la théorie de la fonctionnelle de la densité DFT. 


Jusqu'ici, l'approche de Kohn-Sham est exacte, par comparaison avec 
l'approximation de Born-Oppenheimer. La structure électronique de la 


plupart des matériaux pourrait être résolue très précisément si le terme 
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d’échange-corrélation était connu avec précision. Le problème à ce stade est 
que nous ne connaissons pas encore le terme de d'échange-corrélation, en 
particulier celui qui modélise avec précision les systèmes qui contiennent les 


otbitales atomiques ‘d’?. 
IT.6 Les fonctions d'échange et de corrélation. 


La détermination de la fonction d’échange-corrélation (désignée par XC) est 


la chose la plus importante à savoir pour l’applicabilité de la DET. 


Il existe de nombreuses fonctions différentes, selon les systèmes et les 
propriétés étudiées. Chaque fonction est valable dans des conditions 
particulières et possède ses avantages et des inconvénients ; elle à une forme 


générale qui est la suivante. 


E,cLP(F)] = Î pPex(F)d? (1:17) 


Où 
" E,.[p(r)] est l'énergie totale 
" £,.(r) est la densité d'énergie. 


Différentes fonctions peuvent être caractérisées par la façon dont elles 


échantillonnent la densité entourant chaque électron afin de construire &,.(r). 
Le résultat simple pour la densité d'énergie XC est : 


Rd eue (II.18) 
exc) => el : 


Lorsqu'il est inséré dans (11.17), on obtient : 


Dr 11.19 
Exclp(P)] — 5 ['PGar Nec) je ( ) 


FA 
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Tout simplement, si n,.(”,r”) est connue, l'ensemble du problème de plusieurs 
cofps serait résolu. L'attention est maintenant tournée vers une discussion sur 


les fonctionnalités de XC. 
ITL.6.1 Approximation de densité locale (LDA) 


La fonction la plus élémentaire est l’approximation de la densité locale (LDA) 
proposée par Hohenberg et Kohn La densité d'énergie XC est approximée 
comme celle d'un gaz d'électrons homogène de même densité que le système 


réel, exe [@] [7]: 
L'énergie d’échange-corrélation donnée par la LDA est la suivante : 


Eïc LP(")] = Î p(F)e 0" [o(F)]df (IL.20) 


ex" pour le gaz d'électrons homogène est connue numériquement, c’est une 


hypothèse de départ raisonnable et la LDA fonctionne relativement bien pour 
les systèmes dont la densité varie lentement. La prochaine étape dans 
l'avancement de la fonction énergie XC consistera à incorporer le gradient de 


la densité. 
IL.6.2. Approximation du gradient généralisé (GGA) 


Dans cette approximation GGA, l'énergie d'échange et de corrélation est 
donnée en fonction de la densité électronique et de son gradient. Le gradient 
de la densité peut être incorporé comme un facteur d’amélioration qui 


modifie £XOM, La forme générale peut être écrite comme suit : 


Exe LPGI = Î pe" PIE Lo), vod (20 


En général, les approximations GGA et LDA fonctionnent très bien pour de 
nombreux systèmes. Cependant, elles ne parviennent pas à bien prédire 
certaines propriétés électroniques tel que le gap énergétique qui est sous 


estimé, surtout quand il s’agit de systèmes hautement corrélés. 
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IL.6.3. Potentiel de Becke et Johnson modifié mBJ 


Comme mentionné, la LDA et la GGA sous-estiment le gap énergétique. 
Becke et Johnson [8] ont développé une méthode basée sur la méthode du 
potentiel optimisé (OPE) qui est une méthode très exigeante en calcul mais 


qui améliore considérablement la valeur du gap énergétique. 


La corrélation GGA ou LDA est utilisée conjointement avec l'échange BJ. 
Tran et Blaha [9,10] ont modifié la fonction BJ (MBJ), ce qui a encore amélioré 


les gaps énergétiques ; c’est cette forme qui est utilisée dans cette recherche. 


Le potentiel mB] peut être défini comme : 


ed) = CUER(r) + (3c — 2)- se Pre 
9 (11.22) 


Où : 
Do = De lui l'est la densité électronique. 
 — 1}, pe VE Vis est la densité d'énergie cinétique. 
“ uÆr est le potentiel de Becke-Roussel (BR). 


" cest un paramètre dépendant du système, avec c=1 correspondant au 


potentiel original de Beck-Jonson. 


Pour les matériaux cristallins, Tran et Blaha proposent de déterminer € par la 


relation empirique suivante : 


I7p,. (rl er) (IL.23) 


= a+ pe pan leu GS p,(r') 


Où : 


Voeu est le volume unitaire de la cellule. 
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“ a = —0.012 et B= 1.023 Bohr sont des paramètres ajustés en fonction 


de valeurs expérimentales. 
II.7 La méthode des ondes planes augmentées (APW) 


La solution de l'équation de Schrüdinger avec un potentiel nul est une onde 
plane. Les électrons éloignés des noyaux se comportent comme des particules 
libres et peuvent être représentés par des ondes planes. Les électrons proches 
du noyau (niveaux de cœur) se comportent comme dans un atome libre et, 
par conséquent, des fonctions analogues à celle de l’atome peuvent être 
utilisées pour les représenter. Ce principe est adopté dans la méthode des 


ondes planes augmentées (APW) [11]. 


La méthode (APW) partage l'espace en deux régions, comme indiqué sur la 
figure 11.2. La première est celle des sphères muffin tin sur laquelle la fonction 
d'onde est développée et elle est de type sphérique et la deuxième région est 
appelée espace interstitiel représenté par une fonction d’onde plane. L'onde 
plane augmentée (APW) utilisée dans l'expansion de W£(r, E) est définie selon 


l'équation (IL.24) : 


UE 
/ ÈT B 
| ie. _| T il — , S œ 
| 7 NN 

.  _. Dé j | "+ T | 
|  — | 

\ NRe 

I Le M # 


Figure IL.2 : Méthode FP-LAPW : l’espace est divisé en deux régions : région 


des sphères muffin tin et région interstitielle. 
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Es eiR+R)r 
WV rel 


BEÉCRE)= aliRr : 
D'AEpen Se (29 


im 
im 


e Où 

e kest le vecteur d'onde. 

e Kest le vecteur du réseau réciproque. 

e rest le vecteur de position. 

eV représente le volume de la cellule unité. 


L'ensemble de base APW est k dépendant. 


ASF et E sont des paramètres à déterminer. u} (r',E) sont les solutions de la 
partie radiale de l'équation de Schrôdinger pour un atome libre aux énergies E. 
D'autre part, cela implique que la fonction d'onde à l'intérieur de la sphère de 


muffin tin dépend du paramètre d'énergie E. 


La fonction d’onde ne peut pas être discontinue à la limite de la sphère. La 

continuité de cette fonction d’onde est nécessaire à la frontière séparant la 

sphère du reste de l’espace. Les conditions aux limites de cette frontière 
| , Re de 

permettent de déterminer les coefficients A4%,. L'énergie E est la valeur propre 

de l'équation de Schrôüdinger ; elle est assujettie à une amélioration à l’aide des 


cycles auto cohérents (SC). 
II.8 Méthode d'ondes planes augmentées linéarisées (LAPW) 


L’inconvénient de la méthode APW est la construction des fonctions uf(r',E) 
avec les valeurs propres à rechercher E = ef. Ceci a été résolu par la méthode 


LAPW (par Taylor) en développant uf autour de l'énergie E. 
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Our (r',E) 


ui (r,E) = uf(r',Eo) + (Eo — En) —# 


+ O(E — €i)? 
(1.25) 


L’équation (11.25) permet une grande flexibilité dans le choix de E;. Les deux 
premiers termes sont les paramètres supplémentaires inclus dans l’extension 


APW pour former un ensemble de la base LAPW à une énergie E, fixe. 


1 > 2 
eitE+K)T 


WE (RE) = Wv re 
K d : k K al 1 s k K A ; ” > 
> fa Her, E) + BORN Cr EVE ) TES (IL.26) 
im 
Où 
" (Eo— Ex) = ei et 
du?(r',E ’ ; 
: Le = ü (r", Eo). 
Les coefficients AE et Be sont déterminés à partir des conditions aux 


limites. Les principaux états électroniques à haute énergie et proches des états 
de valence sont appelés états semi-centraux. Le problème avec ces états est 
qu'ils sont fortement liés et vont parfois au-delà de la région des sphères 
muffin tin. Les états ayant le même nombre l, mais différent du nombre 
quantique principal (états de valence et semi-noyau) doivent être traités avec 
précaution en utilisant la méthode LAPW + LO où LO signifie une orbitale 
locale de valeur nulle dans la région interstitielle et dans la région des sphères 
muffin tin ; en d’autres termes, c’est une orbitale localisée dans la sphère 


muffin tin. 


La mise en œuvre de cet ensemble se fait en incluant un terme radial 


j : (R+K ne s ; 
supplémentaire CA ur (r',Ejn dans la région des sphères muffin tin. Tout 
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compte fait, nous avons trois termes à déterminer en utilisant les conditions 


aux limites et la normalisation. 


Semblable à la LAPW + LO), la base des fonctions d’ondes de la APW à été 
améliorée en incluant le concept des énergies fixes et des orbitales locales, 
conduisant à une méthode baptisée APW + lo où le terme lo est ajouté pour 
renforcer la base des fonctions d’ondes de la APW mais là aussi la première 


dérivée est discontinue comme dans la APW. 


II. Méthode des Ondes Planes Augmentées Linéarisées à Potentiel 


Total [FP L/APW/] 


La méthode FP-LAPW [12] combine l'ensemble de la base LAPW avec un 
potentiel complet et une densité de charge du système. Cela signifie qu'on ne 
recouft pas aux approximations qui prennent en compte le potentiel 
interstitiel (V,°) et le potentiel muffin tin (V$-(7)) en raison de l'inclusion d'un 
potentiel enveloppé YxVÉ e;,: et les termes non sphériques dans la région des 
sphères muffin tin comme : 


D Ven au (P) 


L,M TES y 


D Vie Tel 
K 


V(F) = 


(IL.27) 


II.10 Structure du code Wien2k 


Tous les calculs sont effectués principalement avec le code WIEN2kK [13]. 
Dans cette section, nous allons brièvement présenter ce programme. Le code 
WIENZ2Kk est basé sur la méthode à potentiel complet (L) APW + lo pour 


décrire des matériaux cristallins ; il est divisé en deux parties : 
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IL.10.1 Initialisation 


L'initialisation comprend la configuration du fichier de la structure cristalline 
et la vérification de ses informations relatives à la symétrie avec SGROUP et 
SYMMETRY. A l'étape qui suit, à l'aide de LSTART, les densités atomiques 
sont générées à partir des fichiers d’entrée. La densité de charge de départ 
pour le cycle SCF est générée sous la forme d'une superposition de la densité 


atomique à l'aide de la routine DSTART. 
IL.10.2 Calcul auto cohérent 


Dans le calcul auto cohérent, la routine LAPWO génère tout d’abord le 
potentiel de calcul à partir de la densité de charge, d’autre part la routine 
LAPW1 calcule les vecteurs propres et les valeurs propres de la bande de 
valence en diagonalisant la matrice. Habituellement, c’est la partie la plus 
fastidieuse dans le calcul. Dans l'étape qui suit, la routine LAPW2 génère la 
densité de valence à partir des vecteurs propres, LCORE calcule les états et la 
densité du noyau pour un potentiel de symétrie sphérique, puis le mélangeur 
mélange les densités d'entrée et de sortie pour générer l'entrée du cycle suivant 


et contrôle les critères de la convergence. 
II.11 Thermodynamique atomistique ab-initio. 


La DET est connue comme une technique qui traite des systèmes dans leur 
état fondamental, où la température est prise égale à zéro. Souvent les études 
faites sur des matériaux doivent être réalisées en considérant des conditions 
normales de température et de pression. En thermodynamique, les paramètres 
les plus courants sont la température (T) et la pression (P) et toutes les autres 
grandeurs physiques y dépendent. Pour étendre le pouvoir prédictif des 
calculs ab initio aux températures et pressions, le formalise ab initio de la 
thermodynamique à été formulé. L'idée de base est d'interpréter les valeurs de 
la DFT en termes de potentiel thermodynamique. Si l’un quelconque des 


potentielles thermodynamiques est connu à partir des calculs DFT, dans le 
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cadre existant de la thermodynamique, d’autres propriétés peuvent être aussi 


calculées. 


Dans cette partie nous allons présenter les techniques développées pour le 
traitement des propriétés thermodynamiques des systèmes solides selon la 


nouvelle version du modèle quasi-harmonique du code Gibbs 2. 
IL.11.1 Introduction au programme Gibbs2 


Gibbs2 est un programme de calcul qui prend en compte la dépendance en 
pression et en température des propriétés thermodynamiques en phase solide 
à partir de données ab initio, dans le cadre de l'approximation quasi 
harmonique. Le prédécesseur de ce code est Gibbs dont les auteurs sont M. 


À. Bianco, E. Francisco et V. Luaña [14]. 


Dans un calcul typique, pour une phase donnée, l'utilisateur sélectionne une 
maille de volume englobant la géométrie à l’équilibre. À ces volumes fixes, les 
autres paramètres structurels sont relâchés et une courbe E (V) est obtenue, 
c’est l'énergie statique. Dans Putilisation la plus simple possible de gibbs2, 
c'est la seule donnée nécessaire pour générer les propriétés 
thermodynamiques à pressions et températures arbitraires. Dans les cas les 
plus complexes, la structure de bande vibratoire complète de chaque volume 
est lue et les expressions exactes des formules quasi-harmoniques sont 
utilisées. 

Dans un contexte théorique, l'effet de la pression est expliqué de manière 
simple en ajoutant un terme +pV. Cependant, l'effet de la température 
nécessite un modèle thermique : une manière d’inclure la contribution 
thermique des degrés de liberté des cristaux à énergie libre. Une fois le 
modèle thermique choisi, il est possible de construire l'énergie libre de 
Helmholtz vibrationnelle à n'importe quelle température (Fy:2) et de trouver 
ainsi le volume d’équilibre à température (T) et pression (P) données en 
minimisant l’énergie libre de Gibbs hofs d’équilibre qui est donnée par : 
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G*(V,P,T) = E(V) + PV + Fy(@n(V):T) (11.28) 


=  E(V) est l'énergie totale par cellule du composé. 
" D, est la température de Debye. 
"  Fy5 est le terme vibratoire de Helmholtz, qui peut être écrit comme : 


90) 


8T 


(IL.29) 


—0 
Avyip (@»; T) = nKkT | + 3In(1 + e Dr es D(85/T) 


Où 
" nest le nombre d'atomes 
" _ D(8,/T) représente l'intégrale de Debye, 


Le volume d'équilibre V (P,T) est ensuite utilisé pour calculer le reste des 
propriétés thermodynamiques. Dans cette formulation, les degrés de liberté 
internes (c'est-à-dire ceux qui déterminent la géométrie en dehors du volume) 


sont supposés demeurer inchangés du fait des effets de vibration. 


Par conséquent, nous pouvons trouver la courbe du volume d'équilibre 
V (P,T) en minimisant la fonction G*(V,P,T) par rapport au volume V : 


(CC P, 2) . (IL.30) 
ov PT 


De plus, en résolvant l’équation (11.30), on peut obtenir le module de 
compressibilité isotherme BT(T,P), les capacités thermiques Cy et C?, le 


coefficient de dilatation thermique &y [13,14] et l'entropie S°: 


0G*(V;P,T 
By(T,P) = v( ( ). (1.31) 


oV 
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: 0\_30/T (IL.32) 
Cy = 3nky 140 (>) ES + 
Cp = Cy(l + ayYenT) (IL.33) 
0 
S=nk 4 (>) — 3in(1— er) CES 
_ Ven Cv (1.35) 


a 
ur: 


Les grandeurs thermodynamiques ont été étudiées dans le chapitre IV, pour 


prédire les propriétés thermodynamiques de nos composés DMS. 
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CHAPITRE III. PREDICTION DES PROPRIETES 
PHYSIQUES DU COMPOSE « MgSe » DOPE AU 
MANGANESE (Mn) DANS LA PHASE ZINC-BLENDE. 


III.1 Introduction 


La conception et la fabrication de nouveaux matériaux aux propriétés 
intéressantes est un domaine de recherche très actif et très intense. Sur la base 
de ce qui a été expliqué dans le premier chapitre, les semi-conducteurs dopés 
avec un élément magnétique (DMS) représentent une classe importante de 
matériaux incontournable dans l'industrie des nouvelles technologies. La 
recherche sur ces matériaux s’oriente souvent dans le sens de l'exploration et 
du contrôle de leurs propriétés électroniques et magnétiques en fonction de 
la quantité de lPélément magnétique ajouté dans le semi-conducteur 


constituant la matrice hôte. 


Beaucoup de recherches ont été réalisées sur des matériaux DMS cristallisant 
dans la structure zinc-blende avec comme dopant l'élément de transition Mn 
[1-5]. Dans ce type de matériaux, le comportement dépend de la polarisation 
du spin ; en effet, il a été constaté qu’ils se comportent comme des semi- 
conducteurs pour les deux directions de spin, faisant de cette impureté (Mn) 
un bon choix pour avoir une excellente propriété magnétique, ce qui justifie 
notre choix de Mn comme dopant dans l'étude du semi-conducteur binaire 


Mgse. 


Bien que le MoSe et ses alliages ont été abondamment étudiés et que la plupart 
de leurs propriétés sont considérées comme bien connues [6], il demeure que 
ce composé dopé au manganèse (Mn) n’a pas fait objet d’études 
expérimentales ni théoriques, ce qui nous a motivé à prendre en charge son 


étude au plan théorique. 


59 


CHAPITRE III. PREDICTION DES PROPRIETES PHYSIQUES DU COMPOSE 
«MgSe » DOPE AU MANGANESE (Mn) DANS LA PHASE ZINC-BLENDE 


Compte tenu de l'importance des composés Mai.Mn,Se dans diverses 
applications magnétiques et optoélectroniques, nous avons fixé comme 
objectif l’étude théorique de leurs propriétés structurales, électroniques et 
magnétiques dans la structure ferromagnétique zinc blende (B3) via des 


calculs ab initio. 
Les résultats de ce chapitre ont été publiés dans la référence [7]. 
III.2 Considérations techniques 


Les calculs sont réalisés en se basant sur la théorie de la fonctionnelle de la 
densité en spin polarisé SP-DFT [8]. Nous avons utilisé la méthode des Ondes 
Planes Augmentées Linéarisées à Potentiel Total [FP L/APW] [9], 
implémentée dans le code WIEN2k/[10] pour le calcul des propriétés 
physiques des composés Mai ,Mn,Se dans la phase zinc blende (ZB) dont le 


groupe d'espace est F-43m. 


Pour se faire, nous avons utilisé une super cellule standard (1 X 1 X 1) de 
symétrie cubique contenant huit atomes Mo:Mn;Se. Les composés Mo: 
«Mn,Se sont obtenus par remplacement les atomes de Mo par des atomes Mn 


de façon à réalisé les concentrations suivantes : x = 0.25, 0.50, 0.75. 


L’approximation du gradient généralisée, développée par Wu et Cohen WC- 
GGA [11], utilisée pour estimer le potentiel d'échange-corrélation, sert à 
calculer les paramètres d'état fondamentaux, tel que la constante de réseau &o 


et le module de compressibilité B. 


Pour les propriétés électroniques et magnétiques, c’est l'approche de Becke et 
Johnson modifiée par Tran-Blaha (TB-mBJ) [12], associée à approximation 


PBE-GGA proposée par Perdew, Burck et Ernzenhorf [13] qui est utilisée. 


Dans le cadre de la méthode [FP L/APW/], les fonctions d'onde sont 
développées sur une base d’harmoniques sphériques à l’intérieur des sphères 


Muffin Tin non chevauchées entourant les sites atomiques et sur une base 
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d'ondes planes dans le reste de l’espace (espace interstitiel) de la maille 


unitaire. 


On limite le nombre quantique ofbital ! du développement de la fonction 
d’onde dans la sphère atomique à Lx = 10. Le développement de Fourier 
de la densité de charge est limité au vecteur d’onde du réseau 


réciproque Gnax = 12 (a.u) *. 


Le paramètre de convergence Rx X Kmax qui contrôle la taille de 
l'ensemble de base dans le calcul est pris égale à 8 (Kmyax est la valeur 
maximale du module du vecteur d’onde K du réseau réciproque et R,4 est 
le rayon moyen des sphères muffin-tin et dont les valeurs respectives relatives 


aux atomes Mo, Mn et Se sont 2, 2.16 et 2.06 unités atomiques (a.u)). 


Pour le calcul des propriétés électroniques, l’espace réciproque est 
échantillonné par un maillage de 12 x 12 X 12 suivant la procédure de 
Monkhorst-Pack, équivalent à 2000 k points, il a été utilisé pour l'intégration 
de la zone de Brillouin (B2Z). Le processus d’itération est répété jusqu’à ce que 
l'énergie totale calculée du cristal converge avec une précision de 107* Ry 


lorsque la charge du système est convergée à 0.0001 e. 
IIL.3 Résultats et discussion 
IIL.3.1 Optimisation de la structure 


L'étude des propriétés structurales est une étape fondamentale, car elle nous 
permet d'optimiser les grandeurs nécessaires aux calculs des autres propriétés. 
Comme nous l'avons mentionné précédemment, nous ne disposons pas de 
données expérimentales publiées des composés Ma:,Mn,Se avec x —0, 0.25, 
0.50, 0.75 et 1. Par conséquent, nous avons calculé la constante du réseau de 


nos composés à partir la loi de Végrad obéissant à la formule suivante : 


Amg:_xMnySe (x) = (1 — X)amgse + XAmnse CHE) 
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Le paramètre de réseau du composé se présente alors sous la forme d'une 
combinaison linéaire des constantes de réseaux des deux composés binaires 
et sa valeur calculée est ensuite utilisée comme paramètre d’entrée dans le 


code de calcul WIEN2K. 


Nous rappelons que les deux composés binaires MgSe et MnSe ont la 
structure de zinc blende (B3) et leurs paramètres de mailles respectifs sont 
5.89 (À) [14] et 5.90 (À) [15] ; quant aux atomes Me et Mn, ils occupent les 


positions respectives (0, 0, 0) et (74, 74, /4). 


Ainsi, pour déterminer les paramètres structuraux à l’équilibre (paramètre de 
maille a, et le module de compressibilité B,, on minimise l'énergie totale par 
rapport au volume et on l'ajuste à l’aide de l'équation d'état empirique de 


Birch-Murnaghan [16] dont l'expression est la suivante : 


= Et (ep) (re (1-5)+ (0-1 
de PCR" — 1) ue V V (LIL.2) 
Où 
"  V représente le volume de la maille unitaire à équilibre. 


Le module de compressibilité est déterminé au minimum de la courbe E(V) 
par la relation (IIL.3) ci-dessous : 


02E (IIL.3) 


Bo = V — 
L av? 


Les courbes calculées exprimant la variation de l'énergie totale en fonction du 
volume des composés Mo:,Mn,Se (x 0, 0.25, 0.50, 0.75 et 1) sont illustrées 
sur la figure IIL.1. 
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Figure IIL.1 : Variation de l'énergie totale en fonction du volume des 
composés Moi.Mn,Se (x =0, 0.25, 0.50, 0.75 et 1). Les courbes sont calculées 
à l’aide de lapproximation WC-GGA. 


Les valeurs calculées des paramètres structuraux ainsi que d’autres valeurs 
expérimentales et théoriques relevées de la littérature, sont rangées dans 


le tableau IIL.1. 
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Le tableau III.1. Valeurs calculées de la constante de réseau &@o, et du 
module de compressibilité B, des composés Mai,Mn,Se (x —0, 0.25, 0.50, 
0.75 et 1) calculées par lapproximation WC-GGA. Les valeurs sont 


comparées à d’autres valeurs théoriques et expérimentales de la littérature. 


Constante de réseau ao(À) Module de 
compressibilité B (GPa) 
Composé Notre Autres Valeurs Notre Autres 
travail valeurs expérime travail valeurs 
calculées.  ntales calculées. 
MgSe 504 595117]: 589 MA 47.76 48.13 [17] 
5.92 [18] 49.00 [15] 
5:90"[19] 45.86 [19] 
Mgo75Mno25Se 5.90 - - 49.72 - 
MgosoMnosoSe 5.87 - - 52.17 : 
Mpgpo25MnozsSe 5.86 - - 55.60 - 
MnSe 5.85 5.84[20] 590115] 57.86 42.79 [20] 


130.88 [21] 


Les résultats regroupés dans le tableau IIL.1 montrent que nos valeurs de la 
constante de réseau des binaires MoSe et MnSe sont très proches des valeurs 
expérimentales [14,15], et sont également en bon accord avec les valeurs 


théoriques de la littérature [17-20]. 


Sur la figure IIL2, nous avons représenté la variation du paramètre de maille 
a, en fonction de la concentration x des composés Mo1.,Mn,Se. Les valeurs 
de assont comparées à celles obtenues par la loi de Végard (lignes en 
pointillés). De même, nous avons représenté sur la figure IL.3 la variation du 
module de compressibilité de nos composés en fonction de la teneur du 
manganèse, les valeurs de B, sont dans ce cas aussi comparées à celles données 


par la loi de Végard. 
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5.94 —@- Mg, Mn, Se 
Loi de Végard 

5.92 | 

< 5920! 

œ 
5.88 | 
5.86 | 
0.00 0.25 0.50 0.75 1.00 


concentration (x) 


Figure IIL2 : Variation du paramètre de maille a (À) des 
composés Mo: .Mn,Se en fonction de la concentration x. Les calculs sont faits 


en utilisant approximation WC-GGA. 


T é T ÿ T 
—@— Mg, Mn,Se 
: Loi de Végard 
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0.00 0.25 0.50 0.75 1.00 
Concentration (x) 
Figure IIL.3 : Variation du module de compressibilité B, (GPa) des composés 
Mg::Mn,Se en fonction de la concentration x. Les calculs sont faits en 


utilisant l’approximation WC-GGA. 


Sur la figure IIL.2, on observée un écart entre la courbe obtenue par la loi de 
Végard et celle obtenue à partir des valeurs calculées sur la base du paramètre 
du réseau figurant dans le tableau IIL.1. L'écart est plus important pour 
x=0.50, il est relativement faible pour les autres concentrations. On voit 
clairement une diminution dans la valeur de la constante du réseau quand la 


concentration de Mn augmente dans les composés Mai,Mn,Se (x — 0.25, 
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0.50, et 0.75) à cause de la valeur du rayon ionique de Mn qui est plus petite, 


comparée à celle de l'atome Mg. 


Contrairement à ce qui a été observé sur le paramètre du réseau, la valeur du 
module de compressibilité, augmente dans les composés quand la 
concentration de Mn augmente (voir la figure IIL.3). L'écart relevé entre les 
valeurs fournies par la loi de Végard et celles que nous avons calculées est 


faible. 
IIL.3.2 Propriétés électroniques 


Certaines propriétés physiques des solides sont directement liées à leur 
structure de bande d’énergie électronique et à leur densité d'états électronique. 
La connaissance de ces deux propriétés est nécessaire à lidentification de 
matériaux qui peuvent servir efficacement dans les dispositifs 


optoélectroniques, magnéto-optiques et électromagnétiques. 


Dans cette section nous présenterons les structures de bandes d’énergie 
électroniques ainsi que les densités d’états électroniques totales et partielles de 
nos composés, calculées en utilisant l’approximation TB-mBJ, combinée avec 


Papproximation PBE-GGA. 


Il convient de rappeler que la combinaison de approximation TB-mBJ] avec 
lapproximation du gradient généralisé (GGA) ou avec l'approximation de 
densité locale (LDA) s'est avérée efficace notamment dans le calcul des gaps 
énergétiques des matériaux à base de semi-conducteurs [12] dopés aux métaux 


de transition [22-26]. 
IIL.3.2.1 Structures de bandes d’énergie électroniques 


Pour prédire les coufbes de structures de bandes d’énergie électroniques des 
composés Ma,Mn,Se (x — 0 ; 0.25, 0.50, 0.75 et 1) dans la structure zinc 
blende ferromagnétique, nous avons utilisé les valeurs calculées des 


paramètres structuraux discutées précédemment. 
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Nous rappelons que pour le cas du composé binaire MoSe (x = 0), nous avons 


calculé et tracé la structure de bande d’énergie dans la phase non magnétique 


Les résultats obtenus relatifs aux composés Moi,Mn,Se (x = 0, 0.25, 0.50, 
0.75 et 1) sont illustrés dans les figures IIL.4, IIL5, IILG, IIL.7, et IIL8. Le 


niveau de Fermi (E}) est fixé à l'énergie zéro pour servir de référence. 


TB-mBJ 
14 1 
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Figure IIL.4 : Structure de bande d’énergie du binaire MoSe (x — 0) dans la 


phase non magnétique. 
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Energie (eV) 


Figure IIL5 : Structures de bandes d’énergie polarisées en spin [spin 


majoritaire (up) et spin minoritaire (dn)] du composé Mgo7:5Mno5Se. 


Energie (eV) 


Figure IIL6 : Structures de bandes d’énergie polarisées en spin [spin 


majoritaire (up) et spin minoritaire (dn)] du composé Mgo.50Mn0.508e. 
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Figure IIL.7 : Structures de bandes d’énergie polarisées en spin [spin 


majoritaire (up) et spin minoritaire (dn)] du composé Mgo2:5MnosSe. 


Energie (eV) 
© à 8 D © N° À Oo 


WE D OX WKW LT X Wk 
Figure IIL8 : Structures de bandes d’énergie polarisées en spin [spin (up) et 
spin (dn)] du composé MnSe (x = 1). 


Les résultats montrent que le composé MgeSe (x — 0) dans la phase non- 
magnétique (figure IIL.4) possède un gap énergétique direct dans la direction 


T'- T de haute symétrie dans la première zone de Brillouin. 
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Pour les composés de concentrations x = 0.25, 0.50, 0.75 et 1.0, on distingue 


deux cas : 


Cas des spins majoritaires (Up) : dans ce cas, les courbes obtenues sont 
qualitativement semblables à cause de la structure atomique de type ZB 
des composés et à cause de l’environnement chimique des espèces 
atomiques qui ne varie pas significativement. Pour ces composés, on 
observe dans la direction l — lun pseudo gap car, en effet, il est 
difficile d'affirmer l'absence d'états au sommet de la bande de valence. 


La bande d'énergie dans la bande de valence frôle le niveau de Fermi. 


Cas des spins minoritaires (Down) : dans ce cas aussi, les courbes se 
ressemblent pour les mêmes raisons sus évoquées. Nous notons un gap 
direct dans la direction l — F car, dans ce cas, il est bien clair qu'il 
n'existe pas d'états électroniques au niveau de Fermi. Ce résultat 
indique bel et bien le caractère ferromagnétique semi-conducteur de 


nos composés. 


Les valeurs du gap énergétique des différents composés calculés pour les deux 


directions du spin (Up et Down) sont groupées dans le tableau IIL. 


Tableau III.2 : Valeurs du gap énergétique (cas des spins majoritaires ‘up’ et 


minoritaires ‘down’) des composés Mei,Mn,Se (x = 0.25, 0.50, 0.75 et 1.0) 


Gap énergétique (eV) 
EL-T 


Composé En Exp. Autres 
Ce travail (Spin Ce travail (Spin résultats 
Up) Down) 
MgSe 4.12 (NM 4.12 (NM 4.00 [27] 4.10 [19] 
4.20 [29] 
2.60 [19] 
Mgo75 Mno25Se 3.290 3.83 - - 
Mg.50 Mhno50Se 2.682 3:87 = - 
Mpgo25 MnosSe 20) SRE 5 L 
MnSe 1:767 70 2.90 [28] 3:62; [21] 
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Afin de mieux illustrer l’évolution du gap énergétique dans nos composés, 
nous avons tracée sur la figure IIL.9 la courbe donnant sa valeur en fonction 
de la concentration du manganèse (Mn) pour les spins Up et pour les spins 
Down. D'une manière générale, la valeur du gap énergétique décroit avec la 
teneur en manganèse dans les composés mais une différence frappante est 
notée entre les valeurs de cette grandeur quand on change l'orientation du 
spin. En effet, la différence réside dans le fait que la décroissance est beaucoup 
plus importante pour la direction Up que pour la direction down et on peut 
même dire que le gap pour la direction down garde pratiquement la même 


valeur qui est autour de 3.8 (eV). 


4.5 


4.0 Mg, ;.Mn, ,,$e Mg, Mn, .,.Se 
M M S 
_ 3.5 do.50 MN 50? € 
® 
S 3.0 M9, ;Mn, 8e 
LL À 


2:5 


M9 MN 5755 


2-0 F-@- Spin Up 
—Q- Spin Down 


0.00 0.25 0.50 0.75 1.00 


concentration x 


Figure IIL.9 : Gap énergétique des composés Mai.Mn,Se en fonction de la 


concentration de Mn pour les spins up et down. 


Il est possible d'établir une relation donnant la valeur du gap énergétique 
(bande interdite) en fonction de la valeur de la concentration de l’impureté 
dans les composés Mai,Mn,Se (toutes les valeurs de x comprises entre 0 et 
1). Sachant que expression du gap énergétique des composés Ma:,Mn,Se est 


donnée par l'équation (IIL.4) ci-dessous, il est toujours possible d’ajuster cette 
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équation avec les valeurs du gap de la figure IIT-9 pour extrapoler la valeur du 
paramètre (b). 

Egg: _xmnsse)(X) = XEy(mnse) + (À — X)Eg(mgse) — X(1 — x)b (IIL.4) 
Où : 

ME (mg mnxse)(X) est le gap énergétique du ternaire Ma: Mn,Se. 

" E(mgse) Est le gap énergétique du binaire MoSe. 

D: 


g(Mnse) St le gap énergétique du binaire MnSe. 


" best le paramètre de courbure du gap énergétique des composés 


Me: .Mn,Se. 


Les expressions du gap énergétique dépendante uniquement de la teneur du 
manganèse, sont données par les équations (IIL.5) et (IIL.6), ci-dessous pour 


les deux directions de spin up et down respectivement : 


Ep nppe = 4.094 — 3.281x + 0.986x? (LS) 
El eouBp = 4.105 — 1.250x + 0.956x7 (III.6) 


Les équations ci-dessus montrent une faible valeur du paramètre de courbure 
pour les deux directions de spin, up et down : bjŸ = 0.986 (eV) et 
DS = 0.956 (eV). Les gaps énergétiques sont des fonctions décroissantes 


avec la variable x. 
IIL.3.2.2 La densité d’états électronique (DOS) 


La distribution des états électroniques dans les bandes d’énergie d’un solide 
peut être mieux comprise à travers le calcul des courbes de densités d'états 
électroniques (DOS). A cet effet, nous avons calculé les courbes de densité 


d’états électroniques des différents composés et nous les avons représentées 
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sur les figures IIL.10, IL.11, IIL.12, 11.13 et 11.14, dans l'intervalle d'énergie 
allant de -15 (eV) à +15 (eV). 


Le niveau de Fermi EF est mis à zéro et il est indiqué par la ligne pointillée 


bleue verticale. 
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Figure IIT.10 : Courbes de densités d’états électroniques totale et partielle des 
états (2p, 3s) de Mo et (4s, 4p) de Se pour le composé binaire MgSe (x = 0) 
dans la phase non magnétique. 


L'analyse des courbes DOS du composé binaire non magnétique MeSe (x — 


0) peut être faite sur plusieurs intervalles d'énergie : 


a) -11.5 (eV) < E < -10.5 (eV) : cette région est principalement dominée par 


les états ‘4s’ de Se. 


b) -3 (eV) < E < Ek : cette région contient les états électroniques de type ‘4p’ 


de Se. 


c) E > 4 (eV) : cette région est dominée par les états électroniques ‘2p’ et 


‘3s’ de Mg et ‘4p’ de Se. 
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Figure IIT.11 : Courbes de densités d’états électroniques totale et partielle 
polarisées en spin des états (2p, 3s) de Mo, (4s, 4p) de Se et (3d, 3d-t2, 3deç) 
de Mn dans le composé Mgo.75 Mno5Se. 
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Figure III.12 : Courbes de densités d’états électroniques totale et partielle 
polarisées en spin des états (2p, 3s) de Mo, (4s, 4p) de Se et (3d, 3d-t2, 3deç) 
de Mn dans le composé Mgo.50 Mno:508e. 
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Figure III.13 : Courbes de densités d’états électroniques totale et partielle 


polarisées en spin des états (2p, 3s) de Mo, (4s, 4p) de Se et (3d, 3d-t2, 3deç) 
de Mn dans le composé Mgo.25 Mno75Se. 
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Figure IITL.14 : Courbes de densités d’états électroniques totale et partielle 
polarisées en spin des états (4s, 4p) de Se, et (3d, 3d-t», 3d4) de Mn du binaire 
MnSe (x = 1). 
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Selon la théorie du champ cristallin [31], les courbes PDOS (Partial Density 
Of States) des composés Mai,Mn,Se (x= 0.25, 0.50, 0.75 et 1) montrent que 
les états ‘3d° (Mn) subissent une dégénérescence imposée par l’environnement 
tétraédrique (La phase Zinc-blende) de l’atome Mn. Le niveau 3d dans les 


composés est, par conséquent, départagé en deux sous niveaux désignés par 


eg et tz (Annexe 1-A). 


Revenons d’abord à l'analyse des courbes des densités d’états des composés 


Ma.Mh,Se relatives aux figures 111.11, [II.12 et 111.13. 


La première remarque qu’on peut faire sur les courbes DOS est leur 
ressemblance avec, bien évidemment, des petites différences observées sur la 
position, la forme et la largeur des pics. Nous pouvons discuter les courbes 
DOS de nos composés ternaires en divisant l'intervalle -15 (eV) à +15 (eV) 


en trois régions comme suit : 


a) -13.5 (eV) < E <-11 (eV) : cette partie de l’intervalle est formée par les 
ofbitales ‘4s’ de l'atome Se et ceci pour les deux orientations de spin 


(Up et Down). 


b) -5.5 (eV) <E < Er : les états électroniques qui peuplent cette région 
proviennent essentiellement de l'hybridation entre les otbitales ‘3d: , 
‘3d-t de Mn et les orbitales ‘4p” Se, ‘3s Mg pour les spins 


majoritaires ; pour les spins minoritaires, c’est les états ‘4p’ et ‘4s”. 


c) 1.2 (eV) <E <15 (eV): dans cette région, deux intervalles d'énergie 


peuvent se distinguer : 


m 1.2 (eV)<E <5 (eV): sur cet intervalle interviennent les 


états ‘3d et ‘3d-t2° de Mn pour le cas des spins minoritaires. 


mn 5 (eV)<E <15 (eV): les états électroniques appartenant à 


cet intervalle d'énergie sont les états ‘4s’ et ‘4p’ de Se. 
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Enfin, pour le cas du binaire MnSe (figure IIL.14), les courbes DOS peuvent 


aussi être discutées sur plusieurs intervalles comme suit : 


a) -13.6 (eV) <E < -12.7 (eV): cette partie est exclusivement dominée 


par les états ‘4s’ Se pour les spins majoritaires et minoritaires. 


b) -5.4 (eV) <E < (Ek): sur cette région, c’est les états électroniques 
engendrés par l’hybridation entre les états ‘3d.,, ‘3d-t2° de Mn et les 
états ‘4p’ de Se pour les spins majoritaires ; pour les spins minoritaires, 


on trouve uniquement les états ‘4p’ de Se. 


c) 0.5 (eV)<E < 5 (eV): pour les spins minoritaires, cette région ne 


contient que les états électroniques ‘3d. et ‘3d-t7 de Mn. 


d) 5 (eV) < E <15 (eV) : cet intervalle contient les états électroniques ‘4p” 
et ‘4s’ de Se, et 3d4, 3d-t24 de Mn et ceci pour les spins majoritaires et 


minoritaires. 


Par ailleurs, on n’a pas observé la présence d’états électroniques de type ‘d’ de 
manganèse (Mn) au niveau de Fermi et ceci pour toutes les compositions 
étudiées. En effet, aucun pic n’a été repéré (pour le cas des spins majoritaires 
et pour le cas des spins minoritaires) sur les courbes DOS, indiquant le 


caractère semi-conducteur ferromagnétique des trois composés ternaires. 
III.4 Propriétés magnétiques et couplage d'échange 


Les valeurs du moment magnétique total et partiel (sur chaque site atomique) 
des composés Mei,Mn,Se (x — 0.25, 0.50, 0.75 et 1) ont été calculées en 


utilisant l’approximation TB-mB] ; elles sont regroupées dans le tableau III.4. 


Pour rappel, le magnétisme dans nos composés provient de l’impureté (Mn), 
dont le niveau 3d subit un éclatement donnant naissance à deux niveaux 


d'(eÿ,t2,) où, selon la règle de Hund, tous les états électroniques pour les 
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spins majoritaires sont occupés et les états électroniques des spins minoritaires 


sont vides ou partiellement occupés. 


Tableau III.3 : Valeurs calculées du moment magnétique total et des 


moments magnétiques partiels en (ug) dans les alliages Ma.Mn,Se 


(= 0.25, 0.50, 0.75 et 1). 


Moment magnétique (Up) 
Composé M'* M M" M 
MgSe - - - - 
Mgo75s Mnoz5Se 5.000 4.383 0.007 0.027 
Mgo5o MnosoSe 5.000 4.402 0.015 0.057 
Mgo25 MnozsSe 5.000 4.424 0.024 0.075 
MnSe 5.000 4.489 - 0.100 


D’après le tableau II. 3, les valeurs du moment magnétique local sur le site 
du manganèse (Mn) dans les composés ternaires sont inférieures à la valeur 
prévue par la règle de Hund et qui est de 5 (43) constituant le moment 
magnétique total dans le composé. La différence est une conséquence de 
l'interaction d'échange p-d. Par ailleurs, une faible contribution du moment 


magnétique local est induite sur les sites non magnétiques de Mg et de Se. 


Conformément à la théorie du champ moyen, les équations de l'hamiltonien 
présentées ci-dessous [30,31], sont utilisées pour le calcul des constantes 
d'échange Na et NB qui interprètent la force des couplages ‘s-d” et ‘p-d’ 


respectivement. 
Hsa = —XNoa(S)o ŒIL.7) 
Ha = —XNoB{S)o (LS) 
Où : 
" No est la concentration volumique des sites cationiques. 


" xest la concentration en atome de Mn. 
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" «a l’intégrale de recouvrement entre les orbitales s et d. 
" B l'intégrale de recouvrement des orbitales p et d. 
= (S)l'aimantation moyenne par ion Mn. 


Ces deux constantes d’échange peuvent être calculées à partir de la structure 
de bande d’énergie et des propriétés magnétiques en considérant les 


interactions de Kondo [31] habituelles. 


28e (IIL.9) 
Noa = xt) 

_ AEy 1.10) 
NoB nn x(S) 


AE. et AË, sont respectivement données par les équations (IIL.11) et (IIL.12) : 
AE: =Ekie=Els (IL.11) 


ÂËy — Eber + Fa (IIL.12) 


Eÿec €t Elec désignent les minimums de la bande de conduction dans 


les directions respectives down et up. 


Epv et Elsv désignent le maximum de la bande de valence dans les 


directions respectives down et up. 


Les valeurs prédites de AËEc, AËy, Noa et No) des trois composés ternaires et 


celles du binaire MnSe sont présentées dans le tableau IIL.4. 
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Tableau IITL.4 : Valeurs calculées de AËE,, AË,, Noa et NB des composés 
Ma.Mn,Se (x= 0.25, 0.50, 0.75 et 1) en utilisant l’approximation TB-mB]. 


Composé AE AËy Not NoB 
Mgo.75 Mno25Se -0.488 -1.030 -0.781 -1.648 
Mgo.50 Mno.s0Se -0.791 -1.771 -0.632 -1.417 
Mgo25 MnorsSe -0.829 -2.360 -0.442 -1.258 

MnSe -0.228 -2.986 -0.091 -1.194 


En valeur absolue, on note une augmentation des grandeurs AEC et AE, et 
une diminution des grandeurs Na et NB avec l'augmentation de Mn dans les 


composés. 


Le signe négatif affectant NB et Na confirme le mécanisme de double 
échange dans nos composés où les interactions “s-d’ et ‘p-d” sont parallèles et 


confèrent aux composés le caractère ferromagnétique (FM). 


Remarquons également que la valeur de NB est plus négative que celle de 
Na, indiquant que lénergie d'échange est due en grande partie aux spins 


minoritaires. 
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CHAPITRE IV. ETUDE DES PROPRIETES 
STRUCTURALES, ELECTRONIQUES, MAGNETIQUES 
ET THERMODYNAMIQUES DES ALLIAGES Mg:,Cr,Se 

DANS LA PHASE ROCK-SALT POUR LES 
APPLICATIONS SPINTRONIQUE. 


IV.1 Introduction 


Comme nous lavons mentionné auparavant, les semi-conducteurs 
magnétiques dilués (DMS) jouent un rôle important dans le domaine de la 
spintronique parce qu’ils peuvent être magnétiques comme ils peuvent 
posséder d’autres natures selon l'orientation du spin (comportement demi- 
métalliques (HM) : semi-conducteur pour une direction du spin et métal pour 
l’autre direction) [1] tout en présentant une polarisation en spin de 100 % au 
niveau de Fermi (E-). Nous rappelons que, l'intérêt des matériaux HM (Half 
Metallic) dans les composés DMS est incontournable concernant 
l'amélioration de performances de beaucoup de dispositifs ; ces matériaux 
peuvent aussi être exploités dans de nombreuses technologies spintroniques 
telles que : la magnétorésistance géante (GMR), magnétorésistance à effet 


tunnel (TMR), mémoire magnétique à accès aléatoire (MRAM) … etc. 


Du fait de leur sollicitation dans des milieux nécessitant des températures et 
des pressions élevées, les DMS ont également été caractérisés du point de vue 
thermique pour mettre en évidence leurs caractéristiques multifonctionnelles 


[2-5]. 


Des résultats théoriques existent dans la littérature [6-13] sur limpact du 
dopage des semi-conducteurs binaires avec des éléments de transition 
magnétiques dans plusieuts phases cristallines. À titre d’exemples, le composé 
MgsSe dans la phase zinc-blende dopé avec des impuretés magnétiques [14- 


18] a fait l’objet de beaucoup d’études ayant pour objectif la recherche de 
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systèmes (matériaux) avec de bonnes caractéristiques spintronique. 
Cependant, il est à rappeler que la structure cristallographique du MgSe est 
multiformes : zinc blende (B3), wurtzite (B4) et rock-salt (B1) [19] et chacune 
de ces structures est différente des deux autres par le paramètre du réseau et 


par les autres constantes structurales. 


En se basant sur notre recherche bibliographique sur le composé MoeSe, nous 
sommes arrivés au résultat qui privilège l’usage de la structure rock-salt par 
rapport autres structures car elle améliore nettement les propriétés du binaire 
MsSe lorsqu'il est dopé avec certaines impuretés [20-21]. Partant de ce fait, 
nous avons choisi d'étudier le MoSe dans la structure rock-salt dopé au 
chrome et nous avons orienté notre travail de recherche dans le sens 


d’investiguer le comportement demi-métallique dans ce type des composés. 


L'objectif est d'étudier les propriétés structurales, électroniques, magnétiques 
et thermodynamiques des composés Moi,Cr,Se (x = 0.25, 0.50 et 0.75) dans 
la phase rock-salt en utilisant le même type de calcul que les précédents tout 
en considérant en plus le modèle de Gibbs quasi-harmonique pour la 


prédiction des propriétés thermodynamiques. 
IV.2 Considérations techniques 


Nos calculs sont basés sur la théorie de la fonctionnelle de la densité en spin 
polarisé (SP-DFT) [22] et sur la méthode FP-LAPW [23], implantée dans le 
code WIEND2K [24] ; c’est dans ce contexte que les solutions de l'équation de 
Kohn-Sham [25] sont cherchées. 
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Pour simuler les propriétés du semi-conducteur MgSe dopé au chrome (Cr), 
nous avons utilisé une super cellule (1 x 1 x 1) à huit atomes de la structure 
rock-salt (RS) de symétrie cubique dont le groupe d’espace est Fm-3m figure 


IV.1, contenant quatre atomes Mo et quatre atomes Se. 


Figure IV.1 : Structures cristallines (a) jusqu’à (d) des composés Mo, Cr,Se 
(x = 0.25, 0.50 et 0.75). Les sphères grises, bleues et jaunes claires, indiquent 


respectivement les atomes de Cr, Mg et Se. 


Les trois composés ternaires Mo1,Cr,Se (x = 0.25, 0.50 et 0.75) qui ont fait 
l’objet de notre étude sont réalisés par remplacement, dans le binaire MoSe, 
des atomes Mg par des atomes Cr. Le calcul des paramètres structuraux des 
composés est exécuté à l’aide de l'approximation du gradient généralisée WC- 
GGA [26], développée par Wu et Cohen pour estimer le terme d’échange- 


corrélation (XC). 


Pour prédire les propriétés électroniques et magnétiques des composés 
Moi.Cr,Se, nous recourons à l'approche du gradient généralisé PBE-GGA 
[27] de Perdew-Burk-Ernzerhof, combinée avec l’approximation TB-mB] 
[28] de Tran-Blaha modified Becke Johnson. Dans notre cas, les fonctions de 
base Rnt X Kmax Sont choisies égales à 8, où R;+ est le plus petit rayon des 
sphères muffin tin et Kyax est le rayon de coupure pour la base de fonctions 
d'onde. La configuration électronique de Mg est 35°, celle de Cr est 3d°4s! et 
celle de Se est 3d!°4s/4p* ; ces électrons sont traités comme des électrons de 


valence. Les valeurs des rayons atomiques (R;,, muffin-tin) des atomes Mo, 
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Se et Cr sont respectivement 2.31, 2.20 et 2.31 unités atomiques (a.u), alors 
que la densité de charge étendue de Fourier est tronquée à Gux = 12 (a.u)”#, 
et le maillage k de 12 x 12 x 12 équivalent à 2000 K points est sélectionné 
pour la zone de Brillouin (BZ). Pour chaque calcul, la convergence de l’énergie 
totale du système est égale à 10*(Ry) lorsque la charge du système converge 


vets 0.0001 e. 


Enfin, les paramètres structuraux des composés Me:,Cr,Se (x = 0.25, 0.50 et 
0.75) calculés dans la phase ferromagnétique en utilisant approximation WC- 
GGA (à partir des courbes E(V) [29]), sont exploités dans le calcul des 
propriétés thermodynamiques via le modèle quasi-harmonique de Debye, 


implémenté dans le programme Gibbs [30]. 
IV.3. Résultats et discussion 
IV.3.1 Optimisation de la structure 


Afin d’obtenir les propriétés structurales de l’état fondamental des composés 
Moi.Cr,Se (x = 0.25, 0.50 et 0.75) dans la phase rock-salt ferromagnétique, 
nous avons appliqué l'approximation WC-GGA [26]. Les valeurs calculées de 
la constante de réseau à l'équilibre a, et du module de compressibilité BÇ sont 
présentées dans le tableau IV.I, elles sont extraites en ajustant la courbe 
calculée de l’énergie totale du système avec la courbe calculée à partir de 


l’équation de Brirch-Murnaghan [29]. 
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Figure IV.2 : Variation de l’énergie totale en fonction du volume des 
composés Mai,Cr,Se (x = 0, 0.25, 0.50 et 0.75) Les courbes sont calculées à 
laide de l’approximation WC-GGA. 


Les valeurs de nos paramètres structuraux sont comparées à d’autres valeurs 


théoriques et expérimentales relevées de la littérature [19, 31-33]. 


Nous avons également étudié la stabilité des composés Moi,Cr,Se 
(x = 0.25, 0.50 et 0.75) dans la structure rock-salt ferromagnétique en estimant 


leurs énergies de formation E; , obtenues via l'équation suivante [33-35] : 


Er = [Erota(Mg1-xCrxSe) — 1 — XEpux (M9) — XEpuux (CT) (V1) 
— Epux(Se)] 


Où 
Erota(Mg1-xCr,Se) représente l'énergie totale du composé. 


Epux M9), Epux (Cr), et Enux(Se) sont les énergies totales par atome de 


Me, Cr, et Se respectivement. 


" 1—x, x représentent les concentrations d'atomes Mo, Cr dans la maille 


élémentaire. 
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L'énergie de cohésion (E) des composés Ma:,Cr,Se (x = 0, 0.25, 0.50 et 0.75) 
est ensuite prédite à partir de la relation (VI.2) ci-dessous, le but est d'évaluer 


la pérennité de la stabilité thermodynamique des composés : 


Ec = (Erota(Mg1-xCrxSe) == XEotom(Mg) RÉ (OT) (VI2) 
= Etom(Se)|] 


Où 


n  EGtom est l'énergie de l'atome libre dans son état fondamental, elle est 
calculée en supposant que chaque atome est isolé dans une cellule 
cubique de constante de réseau égale à 15 (À) afin de négliger les 


interactions avec les atomes voisins. 
Les valeurs calculées de E et de E; sont aussi indiquées dans le tableau IV.1 


Tableau IV.I : Valeurs de la constante de réseau aÿ(À) du module de 
compressibilité B; (GPa) de l’énergie de formationE;, et de l'énergie de 
cohésion E, des composés Mai,Cr,Se (x = 0, 0.25, 0.50 et 0.75) dans la phase 
rock-salt, calculées par approximation WC-GGA. 


Constante de réseau aç(À) Module de E; Ë 
compressibilité (eV/atome)  (eV/atome) 
B, (GPa) 
Composé Notre Cal. Exp Notre Cal. Exp Notre travail Notre 
travail travail travail 
MgSe 5.46  5.46[19] 5.47[32] 64.62 65.0 [19] - 
B51 [31] 61.0 [33] 

MgrsCro2sSe 5.41 à - 67.69 - - -0.657 -2.928 
MgosoCro.soSe 5.36 - - 73.64 - - -0.528 -3.400 
Mo25CrorsSe 5.31 - - 76.46 - - -0.398 -3.870 


La valeur du paramètre de réseau du binaire MgSe, que nous avons trouvée 


est en très bon accord avec celle obtenue expérimentalement [32], elle est 
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également en bon accord avec les valeurs théoriques obtenues par d’autres 
auteurs [19,31] ; elle décroit avec l'augmentation du taux du dopant (Cr) dans 


les composés à cause de la taille de l'atome Cr qui est inférieure à celle de Me. 


Cependant, la valeur du module de compressibilité dans les composés 


augmente avec l'augmentation du pourcentage de Cr. 


Quant aux valeurs des énergies de formation et de cohésion, elles sont 
affectées d’un signe moins (valeurs négatives) comme il est montré dans 
tableau IV.I. Ces résultats indiquent bel et bien que nos composés peuvent 
ctistalliser avec une bonne stabilité pérenne dans la structure rock-salt 


ferromagnétique. 


Nous rappelons que faute de disponibilité, dans la littérature, de données 
expérimentales ou théoriques du paramètre de réseau et du module de 
compressibilité des composés Ma:,CrSe (x = 0, 0.25, 0.50 et 0.75), nous 
avons donc utilisé nos valeurs calculées de la constante de réseau pour calculer 


les autres propriétés. 
IV.3.2 Propriétés électroniques 
IV.3.2.1 Structure de bande d’énergie et densité d'états électronique 


Dans cette section, nous étudions les propriétés électroniques du composé 
binaire MgSe dopé au chrome (Cr) pour réaliser les composés ternaires 
Moi,Cr,Se (x = 0, 0.25, 0.50 et 0.75) dans la phase rock-salt. Ces propriétés 
apportent des informations importantes sur la nature des composés, sur leur 


utilité et sur leur usage dans les dispositifs spintroniques. 


Nous avons calculé les courbes de structures de bande d’énergie et les courbes 
de densité d'états électronique totale (TDOS) et partielle (PDOS) des 
composés Mo:.Cr,Se à l’aide de Papproche TB-mBJ-PBE, dans la première 
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zone de Brillouin (BZ) et ceci pour les deux orientations du spin Up (vers le 


haut (T)) et Down (vers le bas (J)). 


Les résultats sont présentés sur les figures IV.3, IV.4 et IV.5, où lénergie de 


Fermi (Er) est prise égale à zéro. 


Pour le cas des spins majoritaires (up), les trois composés possèdent une 
nature métallique en raison d’existence d'états électroniques au niveau de 


Fermi « pic centré à (Er) ». 


Mg, 75Cr5 258€ 


0.25 


TB-mBJ_PBE 


AN 


À Semi-conductor | | 
E, =2.667 eV 


Energie (eV) 


DOS (Etats/eV) 


Figure IV.3 : Courbes calculées des bandes d’énergie et de densité d’état 
électronique totale, polarisées en spin du composé Mgo75Cro25Se. Les calculs 


sont effectués via l'approche TB-mBJ-PBE. 
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Energie (eV) 


DOS (Etats/eV) 


Figure IV.4 : Courbes calculées des bandes d’énergie et de densité d’état 
électronique totale, polarisées en spin du composé Mgo.0CrosoSe. Les calculs 


sont effectués via l'approche TB-mBJ-PBE. 
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Figure IV.5 : Courbes calculées des bandes d’énergie et de densité d’état 
électronique totale, polarisées en spin du composé Mgo25CrorsSe. Les calculs 


sont effectués via l'approche TB-mBJ-PBE. 
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Pour les cas des spins minoritaires (Down), les trois composés présentent un 
gap situé entre le haut de la bande de valence (T VB) et le bas de la bande de 
conduction (BCB) dans la direction l' de haute symétrie de la zone de Brillouin 
indiquant leur caractère semi-conducteur. Les valeurs du gap relevées sont 
2.667 (eV), 2.044 (eV) et 1.831 (eV) pour MgrsCrozsSe, MoosoCrosoSe et 
Mgp25Cro.75Se, respectivement. Cette décroissance observée sur les valeurs du 
gap énergétique, avec l’augmentation de la teneur de Cr dans les composés, 
est la conséquence du déplacement du bas de la bande de conduction vers le 
niveau de Fermi toujours dans la direction de haute symétrie de la zone de 


Brillouin. 


Cette double nature (semi-conducteur pour les spins Down et métal pour les 
spins Up) des composés leur confère le caractère demi-métallique et permet 


de les classer parmi les matériaux DMS. 


Une autre grandeur caractéristique des DMS est le gap demi-métallique 
Gam(eV) (Half Metallic gap HM). Ce paramètre est un excellent indicateur de 
qualité des DMS et de leurs applications dans les dispositifs spintroniques; il 
est défini pour une orientation donnée (spin up ou down où on a un semi- 
conducteur) comme étant la valeur la plus petite choisie entre celle 
correspondant à l’écart entre le minimum de la bande de conduction et le 
niveau de Fermi et l’écart entre le maximum de la bande de valence et le niveau 


de Fermi. 


Les valeurs calculées de Gym sont énumérées dans le tableau IV2, et l'on peut 
voir qu’elles deviennent de plus en plus petites avec l'accroissement du taux 


de Cr dans les composés. 
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Tableau IV.2 : Valeurs du gap énergétique (gap direct) Ej7"P°"9) (eV) et du 


gap demi-métallique Gym(eV) des composés ternaires Moi ,Cr,Se. 


Composé Et Poe Gum (eV) 
Mgo.75 CrozsSe 2.667 0.902 
Mpo.50 Cro.508e 2.044 0.382 
Mgp25 CrorsSe 1.830 0.168 


Afin de mieux voir la distribution des électrons dans les bandes d’énergie des 
composés Moo75Cro25Se, MgosoCrosoSe et Mgo25CrorsSe, nous avons tracé 
leurs courbes de densité d'états électronique partielle (PDOS) dans un 
intervalle d'énergie compris entre -15 (eV) et + 12 (eV) et nous les avons 
représentées sur les Figures IV.6, IV.7 et IV.8. Nous rappelons que le niveau 


de Fermi (Er) est toujours mis à zéro (ligne verticale rouge en pointillées). 


Avant de commenter les courbes PDOS, il convient de s’assurer tout d’abord 
d’une caractéristique essentielle des DMS et qui est le degré de polarisabilité 
en spins (P) au niveau de Fermi et dont la valeur est donnée par l'expression 
suivante : 


_ N°(Er) — N°(Er) (V3) 


“ne N'(Er) + N\(Er) 


Où N'(Er) et N'(E-) sont respectivement les valeurs des densités d'états 
électroniques au niveau de Fermi (E-) pour les deux orientations du spin : Up 


et Down. 


Sachant que la densité d’états électronique au niveau de Fermi des composés 
Mg:CrSe (pour le cas des spins minoritaires) est N°'(E;)=0, leur 
polarisation en spins est donc égale à 100 %, d’où la satisfaction de la 


caractéristique fondamentale des DMS. 
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Figure IV.6 : Courbes de densité d'états électronique polarisées en spin : 
totale (TDOS) et partielle PDOS) relatives aux états électroniques (2p, 3s) de 
Meg, (ds, 4p) de Se et (3d, 3d-t4, 3d«) de Cr du composé Mor5CrorsSe. Les 
calculs sont effectués en utilisant l'approche TB-mBJ-PBE. 
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Figure IV.7 : Courbes de densité d’états électronique polarisées en spin : 
totale (TDOS) et partielle PDOS) relatives aux états électroniques (2p, 3s) de 
Me, (ds, 4p) de Se et (3d, 3d-4, 3d-t2) de Cr du composé Mg.:0CrosoSe. Les 
calculs sont effectués en utilisant l'approche TB-mBJ-PBE. 
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Figure IV.8 : Courbes de densité d’états électronique polarisées en spin : 
totale (TDOS) et partielle PDOS) relatives aux états électroniques (2p, 3s) de 
Meg, (ds, 4p) de Se et (3d, 3d-t4, 3d«) de Cr du composé Mo25CrorsSe. Les 


calculs sont effectués en utilisant approche TB-mBJ-PBE. 


Pour les trois composés et dans les deux directions du spin, nous avons mis 


en évidence les contributions des états électroniques à la densité d’états. Dans 
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la région située entre -15.5 (eV) et -12.5 (eV) en dessous de Er, seuls les états 
‘4s’ (Se) existent ; entre -7 (eV) et 3.5 (eV), on trouve les états ‘4p’ (Se) et ‘3d° 


(Ch. 


Concernant l'intervalle où l’énergie est supérieure à 3 eV, on à une présence 


des états ‘34°’ du chrome, des états ‘2p” et ‘3s’ de Mg et ‘4s’ et ‘4p” de Se 


Pour rappel, les états ‘34° (Cr) subissent une dégénérescence et se divisent en 
deux sous niveaux désignés par ‘e et ‘t selon la théorie du champ cristallin, 


en relation avec l’environnement octaédrique de l’atome Cr (Annexe 1-B). 


La valeur de la densité d’états électronique au niveau de Fermi pour le cas de 


spin minoritaire (N°(Ek)) dans les composés, augmente en fonction du taux 
de chrome. Nous avons relevé les valeurs suivantes : 0.61 (états/eV) pour 
Movr5Cro2sSe, 1.13 (états/eV) pour MgosoCrosoSe et 1.86 (états/eV) pour 
MaoCrorsSe. Une forte contribution des états ‘34? (Cr) est observée, elle 
domine considérablement les états ‘4p’ (Se) et ‘3d-t3 (Cr). Pour cette 
ofientation du spin, le caractère métallique (spin Up) des composés résulte de 


la forte hybridation entre les états ‘4p’ (Se) et ‘3d” (Cr). 


Les courbes de densité d’états électronique (cas des spins Down) affichent 
clairement une valeur nulle au niveau de fermi (N'(E-) = 0) et ceci pour toutes 


les concentrations considérées, ce qui donne une polarisation P = 1. 


L'énergie d'échange (A;(d) = d, — dr), définie comme étant l'écart d’énergie 
entre les pics des états ‘3d” du chrome (observés sur les courbes DOS pour 
les deux configurations du spin : Up et Down), a également été calculée. Elle 
a été comparée à l'énergie du champ cristallin ionique 
(Ecrys = di2g — deg)dans le but de mettre en évidence un ordre 


ferromagnétique stable. 
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En outre, l'énergie d'interaction ‘sp-d” (A,(pd) = Ej —Ef) entre bandes de 
valence (spin Up et spin Down) a été calculée dans le but d’apporter un plus 


à la compréhension du ferromagnétisme dans nos composés. 


Selon le tableau IV.3, les valeurs calculées de A,(d) sont plus importantes que 
celles de Etyys = di2g — deg. Ces résultats sont un bon indicateur de l’existence 
du ferromagnétisme dans les composés. Le signe négatif affectant la valeur de 
A,(pd) signifie que le potentiel effectif associé aux spins minoritaires est plus 
important que celui associé aux spins majoritaires. Sur la base de ces résultats, 


le comportement HM des composés est confirmé. 


Tableaux IV.3 : Valeurs de l'énergie d'échange A,(d), de l'énergie du champ 
cristallin ionique Etrysta et de l'énergie d’interactions (pd) des composés 


Mars5Cro25Se, Mgo50CrosoSe et Moo25CrorsSe. 


Composé AEcrystat A,(d) A,(pd) 
Mgr: CrozsSe 0.463 3.918 -0.514 
Mgo.50 Cro.soSe 0.490 4,000 -0.706 
Mgo25 Cro7sSe 0.816 4,544 -0.680 


IV.3.2.2 Densité de charge électronique 


Afin d'étudier les liaisons entre les atomes dans nos composés, nous avons 
calculé leurs courbes de densité de charge électroniques en utilisant 
l'approximation TB-mBJ-PBE. Les calculs sont effectués dans un espace 
bidimensionnel 2D et précisément dans le plan (110) en prenant en compte 
les deux configurations de spins (Up et Down). Les courbes sont présentées 


sur la figure IV.9. 
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ÉPTÉETETITEEI LENTILLE L 


Figure IV.9 : Courbes de densité de charge électronique des composés 
Me;Cri,Se (x = 0.25, 0.50 et 0.75) pour les deux configurations de spins (Up 
et Down), calculées à l'aide de approximation TB-mBJ-PBE. 


Visiblement, la figure IV.9 ne montre pas de différences notables entre les 
courbes de densités de charge électroniques de valence totale et ceci pour tous 
les composés étudiés. Le tracé du contour fait ressortir le caractère ionique de 
la liaison Mg-Se en raison de la différence d'électronégativité de Pauling et 
d'une forte liaison covalente Cr-Se, tandis que la répartition de la charge entre 
les atomes identiques est synonyme du caractère covalent partiel. En d'autres 
termes, la charge de valence autour des sites Mg provient des états s, et celle 


autour des atomes de Se, quant à elle, elle provient des états ‘s’ et ‘p’. 
IV.3.3 Propriétés magnétiques et couplage d'échange 


En utilisant l'approximation TB-mBJ-PBE, nous avons calculé le moment 
magnétique total des composés DMS Me:_,Cr,Se (x = 0.25, 0.50 et 0.75) dans 


la phase rock-salt ferromagnétique. Nous avons également calculé le moment 
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magnétique local sur chaque site atomique pour chacun des composés ; les 


résultats obtenus sont résumés dans le tableau IV.4. 


Nous rappelons que, d’après la règle de Hund, la configuration électronique 
des états ‘3d” (Cr) dans les composés Ma.CrSe est d'(t5,eg); par 
conséquent, la bande ‘3d” (Cr) relative aux spins majoritaires est partiellement 
remplie et contient trois électrons alors que celle correspondant aux spins 


minoritaires est vide. 


Tableau IV.4 : Valeurs calculées des moments magnétiques totaux et locaux 
ainsi que des constantes d'échange des composés Mai,Cr,Se 
(x = 0,25, 0,50, 0,75). Les calculs sont effectués via l'approximation TB-mBJ- 
PBE. 


Composé Mu) Mu) ME) Mug) Noa Nb 


Mpgp.75 CrozsSe 4.000 3.654 0.006 -0.026 3.917 -1.028 
Mgp.50 Cro.s0Se 4.000 3.639 0.013 -0.032 1.810  -0.706 
Mpgo25 Cro7sSe 4.000 3.668 0.021 -0.065 1.012 -0.453 


Les valeurs des moments magnétiques locaux du dopant (Cr), pour toutes les 
concentrations étudiées, telles que indiquées dans le tableau IV.4, sont toutes 
inférieures à la valeur 4 (u3), fournie par la règle de Hund. Les écarts observés 
entre ces valeur et la valeur 4 (u3), sont dus principalement à la forte 
interaction (hybridation) entre les états électroniques de type ‘d”’ de Cr et les 
états électroniques de type ‘p’ de Se conduisant à une diminution de la valeur 
de M'°* (ug) avec l’augmentation du taux de Cr dans les composés. En plus 
du moment magnétique propre au chrome dans les composés, on note aussi 
l'existence d’un moment magnétique localisé de faible magnitude sur les sites 


atomiques non magnétiques de Mg et Se. 
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Notons que le signe (—) qui affecte les valeurs du moment magnétique local 
de atome Se signifie que l'interaction entre son moment magnétique de spin 
et celui de Cr est antiferromagnétique, alors que le signe (+) affectant les 
valeurs du moment magnétique de Mc signifie que l’interaction entre son 


moment magnétique de spin et celui de Cr est ferromagnétique. 


Les valeurs des constantes d'échange Na et NB des composés ternaires sont 
calculées à partir des équations (IV.4) et (IV.5) ci-dessous [40,41] ; elles sont 


portées dans le tableau IV.4. 


_ Ac AV.4) 
Na xt) 

_ AEy (V5) 
NB E x(S) 


Où : 
" No est la concentration volumique des sites cationiques. 
" xest la concentration en atome de Cr. 
= (S) l'aimantation moyenne par ion Cr. 


" AFcet AE, sont respectivement données par les équations (IV.6) et 


ŒV.7) : 
AE = Eve EE EE AV.6) 


AEy = Eve ES: Elev AV.7) 


"  Eÿnc et Emsc désignent les minimums de la bande de conduction dans 


les directions respectives down et up. 


"  Ehpv et Env désignent le maximum de la bande de valence dans les 


directions respectives down et up. 
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Pour les trois composés, les valeurs de Na sont positives et celles de NB sont 
négatives indiquant respectivement que l'antiferromagnétisme et le 
ferromagnétisme se produisent, entre la bande de valence et les états ‘3d° de 


Cr et entre la bande de conduction et ces mêmes états ‘3d’. 
IV.3.4 Propriétés thermodynamiques 


Afin de prédire les propriétés thermodynamiques des composés Mo:-,Cr,Se 
(& = 0,25, 0,50 et 0,75) dans la phase ferromagnétique rock-salt, nous avons 
appliqué la théorie du cristal quasi-harmonique de Debye telle que mise en 


œuvre dans le programme Gibbs [42,43] combinée à la méthode FP-LPW. 


L'étude des grandeurs thermodynamiques macroscopiques de nos composés 
est faite en fonction de la température dans l'intervalle allant de O0 à 1200 


(K°) avec une pression variant entre 0 et 20 (GPa). 
IV.3.4.1 Capacité calorifique 


Les figures VI.10 et VI.11 montrent l’évolution des courbes des capacités 
calorifiques à pression constante C} et à volume constant Cy des composés 


Mai.Cr,Se en fonction de la température et de la pression. 


Les valeurs de ces deux grandeurs thermodynamiques (C, et Cy), enregistrées 
entre 0 et 300 (K), augmentent d’une manière exponentielle selon la loi T* et 


ont des comportements similaires pour les trois composés. 
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Figure IV.10 : Variation de la capacité calorifique à pression constante C} en 
fonction de la température des composés Mo, Cri,Se à différentes valeurs de 


la pression. 
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Figure IV.11 : Variation de la capacité calorifique à volume constant Cy en 


fonction de la température à différentes valeurs de la pression, des composés 


Mpg:CriSe. 


Pour T > 300 K, Cy tend vers une valeur limite connue sous le nom de limite 
de Dulong-Petit ; par contre, sur ce même intervalle de températures, C» ne 


tend pas à une valeur constante. 


Les valeurs de Cy dans les conditions normales de température et de pression 
des composés MosCro25Se, MeosoCrosoSe et  Moo25Cror5Se sont 


respectivement 69.0025 ( mol=!K71), 69.34634 (J mol”"!K71), et 72.1579 
{ mol- KT. 
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IV.3.4.2 L’entropie S 


Les courbes d’entropie S, enregistrées en fonction de la température à 
différentes valeurs de la pression des composés Mg:-,Cr,Se, sont tracées sur 


la figure IV.12. 
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Figure IV.12 : Variation de l'entropie S des composés Ms;Cri,Se en fonction 


de la température, à différentes valeurs de la pression. 


Les courbes montrent une importante augmentation de l’entropie S avec la 
température, mais elle diminue avec l'augmentation de la pression pour une 


température donnée ; ce comportement est observé dans les trois composés. 
IV.3.4.3 Le volume V 


La figure IV.13 présente la variation du volume des composés DMS 
Me;Cri,Se en fonction de la température pour différentes valeurs de la 


pression. 
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Figure IV.13 : Variation du volume V en fonction de la température pour 


différentes valeurs de la pression des composés Me,Cri,Se. 


Les résultats montrent pour une pression P donnée, que V augmente d’une 
manière monotone avec la température, en particulier pour les grandes valeurs 
de P. D'autre part, V diminue avec l’accroissement de P et du taux de Chrome 


(Cr) dans les composés. 
IV.3.4.4 Le coefficient de dilatation 


La variation du coefficient de dilatation thermique a des DMS Ms: Cri,Se, en 
fonction de la température dans l'intervalle borné entre 0 et 1200 (K) a 
également été étudiée et à différentes valeurs de la pression (P — 0, 4, 8, 12, 


14 et 20 GPa). Les courbes de cette grandeur sont tracées sur la figure IV.14. 
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Figure IV.14 : Variation du coefficient de dilatation thermique a en fonction 
de la température à différentes valeurs de la pression pour les composés 


Mpg:CriSe. 
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Pour une pression donnée, « augmente rapidement avec la température dans 
l'intervalle allant de 0 à 300 (K) ; au-delà de 300 (K), il continu à augmenter 
mais lentement et tend vers une valeur presque constante aux grandes valeurs 


de la température. 


Quand la valeur de la pression augmente, on note une diminution de la valeur 
de « pour une composition donnée. Enfin, pour une pression et une 
température données, a décroit avec l’augmentation du taux de Chrome (Cr) 


dans les composés. 
IV.3.4.5 Le module de compressibilité B 


La figure IV.15 montre l'évolution du module de compressibilité B des 


composés Mo,Cri,Se avec la température à différentes valeurs de la pression. 
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Figure IV.15 : Variation du module de compressibilité B des composés 


Me;Cri,Se en fonction de la température à différentes valeurs de la pression. 


Nous observons qu’à une pression donnée, la valeur de B diminue avec la 
température ; par contre quand la pression augmente, B augmente aussi. 
Notons que l'augmentation de la concentration dans les composés n'a pas 


d'effet significatif sur les valeurs de B. 
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L'exploitation du spin de l’électron en plus de sa charge électrique à donné 
naissance à de nouveaux dispositifs dans le domaine du transport et du 


traitement de l’information, c’est le domaine de la spintronique. 


Les semi-conducteurs magnétiques dilués (DMS) occupent une importante 
place dans ce domaine ; parmi ceux-ci, nous citons les semi-conducteurs de 
type II-VI qui sont des composés prometteurs en raison des avantages qu'ils 
présentent notamment leur super efficacité en polarisation et en injection de 
spin ainsi que leur structure cristalline qui est compatible avec les appareils 


existant. 


L'étude menée dans le cadre de cette thèse est théorique et repose sur la 
théorie de la fonctionnelle de la densité polarisée en spin (SP-DET) qui est 
une méthode de calcul ab initio dotée de performantes capacités pour décrire 
efficacement les systèmes magnétiques, en particulier les semi-conducteurs 


magnétiques dilués (DMS). 


Dans ce contexte, nous nous intéressons au calcul des propriétés physiques 
du composé binaire MoSe dopé aux métaux de transition (Mn et Cr) à 
différentes concentrations. Les différentes propriétés sont calculées en 
utilisant la méthode des ondes planes augmentées avec un potentiel total (FP- 


LAPW), implémentée dans le code WIEN2k. 


Dans un premier temps, nous avons étudié les propriétés structurales, 
électroniques et magnétiques des composés Mo:,Mn,Se dans la phase zinc- 
blende en variant la concentration du dopant (x = 0, 0.25, 0.50, 0.75 et 1). 


Nous sommes arrivés aux résultats suivants : 


Les valeurs calculées du paramètre du réseau et du module de compressibilité 


à l'équilibre des composés binaires MoSe et MnSe sont en bon accord avec 
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les valeurs théoriques et expérimentales relevées de la littérature. Dans les 
composés ternaires, il a été observé une diminution de la valeur du paramètre 
du réseau avec l’augmentation du taux de limpureté (Mn) ; inversement, la 
valeur du module de compressibilité augmente quand la concentration de Mn 


augmente. 


L'analyse des bandes d'énergies électroniques et des densités d'états 
électroniques calculées en spins polarisés a montré que les trois composés 
ternaires sont des semi-conducteurs ferromagnétiques à large gap direct dans 
la direction (1 et ceci pour les spins majoritaires (Up) et minoritaires 
(Down). Cependant, les valeurs du gap sont assez différentes pour les deux 
orientations du spin ; en effet, elles varient en diminuant avec l’augmentation 
de limpureté pour les spins majoritaires mais elles sont presque constantes et 


relativement élevées pour les cas des spins minoritaires. 


Ce comportement ferromagnétique semi-conducteur en spins polarisés des 
ternaires peut leur conférer d’excellentes propriétés optiques et magnétiques 


à la fois. 


L'étude magnétique de nos composés à montré que les valeurs calculées du 
moment magnétique total étaient égales à 5 (18), provenant en grande partie 
du dopant (Mn). Nous avons également mis en évidence l’existence de petits 
moments magnétiques (moments de faible magnitude) sur les sites non 
magnétiques de Mg et Se, créés par l’hybridation impliquant les états ‘p’ de 
(Se) et ‘d’ de (Mn). 


Dans un deuxième temps, nous avons orienté notre recherche sur le même 
semi-conducteur binaire (MgSe) mais dopé cette fois-ci au chrome (Cr). 
L'étude avait deux objectifs : le premier était l’investigation du caractère demi- 
métallique (HM) dans les composés ternaires Moi.Cr,Se, absent dans les 
composés Mo:.Mn,Se, et le second était l'étude de leurs propriétés 


thermodynamiques non disponibles dans la littérature. L’étude des composés 
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ternaires Mo1.Cr,Se est faite dans la structure rock-salt pour x= 0.25, 0.50 et 


0575: 
Les résultats réalisés peuvent être présentés comme suit : 


L'étude structurale à mis en évidence une diminution de la valeur du 


paramètre du réseau lorsque la concentration de l’impureté (Cr) augmente. 


L'analyse des densités d’états électroniques (DOS) et des bandes d’énergie 
électroniques a montré que les composés Mo:,Cr,Se, dans la structure rock- 
salt, possèdent un caractère demi-métallique, principale caractéristique des 


DMS. 


Les valeurs du moment magnétique calculées sur les différents sites atomiques 
dans les composés Moi,Cr,Se indiquent une forte contribution provenant de 
l'élément magnétique (Cr) ; les autres contributions provenant des sites 
atomiques non magnétiques de Mg et Se sont relativement faibles et leurs 
ofigines sont expliquées par l'hybridation impliquant les états électroniques de 


type ‘p’ (Se) et de type ‘d’ (Cr). 


Les résultats thermodynamiques fournis par le modèle quasi-harmonique de 
Debye, mis en œuvre dans le code de Gibbs, ont montré que la concentration 
du dopant (Cr) et la pression n’ont pas d’impact significatif sur les grandeurs 
thermodynamiques C, et C, ; par contre, un effet est observé sur les autres 


grandeurs S, V et «. 


Cette étude théorique du semi-conducteur MgSe dopé au manganèse (Mn) et 
au chrome (Cr) est révélatrice, elle à montré que selon la nature du dopage, 
on peut obtenir soit un ferromagnétique semi-conducteur soit un demi-métal. 
Ces résultats sont un bon indicateur aux expérimentateurs en matière 


d'orientation des travaux de recherche sur ces matériaux. 
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Ce travail ouvre des perspectives avenirs étant donné que certains aspects 
n’ont pas encore été abordés. Par conséquent, nous envisageons compléter 


notre étude sur ces matériaux par les travaux suivants : 
Etude des propriétés optiques et élastiques. 


Etude du MgSe dopé, sous forme de couches minces (Etude à deux 


dimensions). 
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Annexe 1 


A. Théorie du champ cristallin dans le cas de la symétrie 


tétraédrique 


Un tétraèdre peut toujours être représenté dans un cube (figure A-1). Dans 
les cas des complexes tétraédriques (champ tétraédrique), les ligands sont 
placés dans les positions entre les axes x, y,z (en réalité ils occupent 4 sommets 


du cube, voir figure ci-dessous). 


. ns os mnt me M 
2/5 À 


L At 


\ -3/5 A 
lon libre Es v 


—— mm à 
Symétrie tétraédrique 
Dédoublement des orbitales d dans un champ tétraédrique 


Figure A-1 : Influence du champ cristallin tétraédrique sur les niveaux des 
otbitales d menant à une levée de dégénérescence et à apparition des niveaux 
tri 

La levée de dégénérescence est traduite par l’écart A; (énergie de 
dédoublement des otbitales d pour une symétrie tétraédrique), séparant les 
otbitales e, et t2,. Dans ce cas-là, les otbitales dx? — y? et dz? sont stabilisées 
(dégénérescence 2, symbole e) en présentant un écart d'énergie entre leurs 
niveaux d’une valeur égale à -=A par électron. Par contre, les 
otbitales dxy, dyz, dxz sont déstabilisées (dégénérescence 3, symbole t,), la 


; er ; N 2 2 2 
déstabilisation va correspondre à un écart +=A; par électron. 


Dans le cas du champ tétraédrique, l'énergie de dédoublement est plus petite 
que dans le cas du champ octaédrique. Pour le même ligand et les mêmes 


distances M — L, le rapport des deux valeurs est 4/9 A;. 
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B. Théorie du champ cristallin dans le cas de la symétrie 


octaédrique 


Lorsque les otbitales d sont dans un champ de symétrie octaédrique (figure 
B-1), on observe une levée partielle de dégénérescence, trois orbitales d 
seront stabilisées et les deux autres seront déstabilisées. La stabilisation ou 
déstabilisation des ofbitales atomiques d dépend de leurs formes et 
orientation (qui dépendent du nombre quantique magnétique m,). Les 
otbitales dx? — y?, dz? ont leurs lobes orientés vers les ligands (charges 
négatives) et comme les orbitales présentent une densité électronique 
importante (charges négatives), leur niveau énergétique va monter 
(répulsion). La situation est inversée dans le cas des trois orbitales 
dxy, dyz, dxz parce que leurs lobes sont orientés vers les directions diagonales 


(entre les ligands). 


Dans ce cas, les orbitales t,, abaissent leur énergie, elles seront donc 


plus stables et cette stabilisation est traduite par un écart d’énergie égale à 


2 > ; : . , 
—=A par électron. Par contre, l'énergie des orbitales e, a augmenté et la 


; e : N 2 3 2 
déstabilisation va correspondre à un écart de +74, électron. 


3 2 à 34? 
A — | / 8, 
VS de 
= ms es em mms Y A 
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Figure B-1 : Effet du champ cristallin octaédrique sur les niveaux des 
otbitales d menant à une levée de dégénérescence et à l'apparition des niveaux 


€g et 29. 
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Abstract 

In this paper, we have investigated the structural, electronic, and magnetic properties of magnesium selenium (MgSe) doped 
with transition metal manganese (Mn) impurity in the cubic diluted magnetic semiconductor (DMS) zinc blende structure. 
The compounds which we are interested are as Mg1_,Mn,;Se where x change between 0 and 1 by step 0.25. AII properties 
are studied, using first-principles calculation of density functional theory under the framework of the full-potential linearized 
augmented plane waver (FP-LAPW). In our study, we employed the Wu-Cohen generalized approximation (WC-GGA) to 
optimize the crystal structure, whereas Tran-Blaha modified Becke-Johnson potential (TB-mBJ) as a new functional was 
applied to compute the electronic and magnetic properties in order to get some better degree of precision. The electronic band 
structures and density of state plots reveal ferromagnetic semiconducting behavior in these compounds, and the exchange 
constants Noa and NoB are calculated to validate the effects resulting from exchange splitting process. Moreover, for each 
concentration x, the value of total magnetic moment has been estimated to equal to 5 8. The important magnetic moments 
values obtained in these compounds indicate the potential for their use in spintronic devices. 


Keywords MgMnSe - Magnetic properties - TB-mBJ : DMS alloys - Exchange mechanism 


1 Introduction [1-3]. To move forward in this new science, we had to find 

semiconductor materials having also magnetic properties, 
The design and manufacture of new materials with often which are called diluted magnetic semiconductor (DMS). 
surprising properties 1s a very active field of research and Where the intensive research of these materials was initiated 
modern technology. In this context, semiconductor doped by the successful doping of non-magnetic semiconductors 
with magnetic impurity represents an important class Of by magnetic atoms like the 3d transition metal (TM) (Mn, 
materials used in the industry of new technologies; this Cr, Ni, Fe, Co, Ti, etc.) or by rare earth metals (RE) (Ho, 
is because it is possible to control the various physical  Gd, Nd, Sm, etc.) [4]. Another important feature of DMSs 
properties such as the electronic and magnetic properties of js that the energy band gap and other physical parameters 
the system by adjusting the content of a component of the can be controlled by varying composition of magnetic ions 
alloy. in these materials [5]. 

In parallel to this important technological development, The main thrust behind this research is to develop new 
spintronic is a new generation of microelectronic, which  spintronics devices like spin valves, spin light-emitting 
exploits the spin of charge carriers in the emerging field of diodes, magnetic sensors, logic devices, and ultra-fast 
promising materials for spin-based multifunctional devices optical switches [61]. 

On the other hand, several experimental [7-13] and 
theoretical [14-21] research is fixed on the DMSs-based 
magnetic element doped II-V and II-VI semiconductor to 
predict their magnetic properties for usages in spintronic 
application. As well, numerous groups of researchers 
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Zni-,Mn,S [23], Cdi-,Mn,Te [24], Be;_,Mn,Se, and 
Be;_,Mn,Te [25]. It was found that these alloys exhibit 
a semiconducting nature in both spin channels, making this 
impurity (Mn) the best choice for finding excellent magnetic 
property, the reason for which it was used in our theoretical study. 

In addition to that, the alkaline earth chalcogenides 
(AECSs) (I = Be, Ca, Sr, Ba; VI = O, S, Se, Te) have 
engrossed much interest, due to their large band gaps and 
low dielectrics constants [26-29], making these compounds 
important for a large variety of optoelectronics devices. 

Specifically, magnesium selenium (MgSe) is an attractive 
binary semiconductor of the II-VI family, which can 
excite in several phases [30], such as zinc blende with in 
experimental wide direct band gap of 4 eV at 0 K [31] and 
lattice constant of 5.89 À [32]. Used in high-temperature, 
high-power blue and ultraviolet wavelength optics [32- 
35] in ZnSe-based laser diodes and light-emitting diodes 
(LEDs) as a cladding material [35]. Furthermore, some 
theoretical studies are available in the literature which 
reported the magnetic properties of MgSe semiconductor 
doped with 3d TM using the first-principles calculation 
[36-391]. 

According to the best of our knowledge, there are no 
experimental and theoretical studies of structural, elec- 
tronic, and magnetic properties of Mn-doped MgSe for 
the full composition range (0-1), which restricts the com- 
parison of our results. In the light of importance of 
these compounds in various magnetic and optoelectronic 
applications, our objective in this paper is to investi- 
gate the structural, electronic, and magnetic properties 
of Mgi-;kMn,Se DMS compounds in the ferromagnetic 
zinc blende (B3) structure by using the ab initio cal- 
culation within the spin-polarization density framework 
(SP-DET). 


2 Theory and Method of Computation 


The calculation of the present study is performed in 
the formwork of the SP-DFT [40] formulation. We have 
employed the full-potential linearized augmented plane- 
wave (FP-LAPW+lo) method [41] as implemented in the 
WIEN2K code [42] to solve the Kohen-Sham equation [43] 
in order to investigate physical properties of Mg1_,Mn,Se 
(0 < x < 1) compounds in the zinc-blende phase (space 
group F-43m 216). 

However, we consider a standard eight-atom supercell 
Mgs-hMn,Se) which correspond to (1 x 1 x 1) of 
zinc-blende phase with cubic symmetry; we obtained the 
structures of Mgi_,Mn,Se by replacing the Mg atoms 
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with Mn to get the compounds with 25, 50, and 75% 
concentrations, respectively. 

For x — 0.25 and x = 0.75, the host lattice exhibit a 
cubic structure with space group 43m, while for x = 0.50 
the crystal present a tetragonal structure with space group 
Päm2. 

The generalized gradient approximation function devel- 
oped by Wu and Cohen (WC-GGA) [44] was used as the 
exchange correlation potential to calculate the ground state 
parameter, which provides an accurate account of structural 
properties of solids [45, 46], where for the electronic and 
magnetic properties the Tran-Blaha modified Becke John- 
son (TB-mBJ) [47] approach combined with PBE-GGA are 
used. It is worth reminding that this approach is not valid for 
the computations of the equilibrium structural parameters. 

As well the combination of TB-mBJ with GGA or 
local density approximation (LDA) have proved to be 
relatively better method for computing energy gaps of sp 
semiconductor, wide band gap semiconductor and transition 
metalbased materials [47-52]. 

Also, we are performed using plane wave cutoff of 
Rmt X Kmax = 8, which control the size of matrix for 
energy convergence, the Rm+t is the smallest radii of muffin- 
tin spheres and Kmax is the cutoff for the wave function 
basis. Whereas the Fourier expanded charge density was 
truncated at Gimax — 12 atomic units (au) !, the values of 
the atomic radius (muffin-tin RMT) were chosen equal to 
2, 2.06, and 2.16 for Mg, Se and Mn respectively in atomic 
units. The maximum angular momentum quantum number 
lmax is Kept at 10 and a k-mesh of 12 x 12 x 12 equivalent 
to 2000 & points was used for the Brillouin zone (BZ) 
integration. Finally, the total energy convergence is equal 
to 107* Ry when the charge of the system is converged to 
0.0001 e. 


3 Results and Discussion 
3.1 Structural Properties 


Firstly, by minimizing the total energy with respect to the 
volume and fitting it to the empirical Birch-Murnaghan [53] 
equation of state, the equilibrium lattice constant a and the 
bulk modulus (B) for Mg1-,Mn,Se alloys are evaluated in 
zinc-blende phase (B3) for x vary between 0 and 1 by step 
of 0.25. 

Both MgSe and MnSe have zinc-blend (B3) structure, 
where Me and Mn atoms are located at (0, 0, O0) and the Se 
atom at (1/4, 1/4, 1/4) position. In order to gain deep insight 
into the physical properties, we have performed volume 
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Table 1 Calculations of the 
lattice constant ao and the bulk 


Lattice constant ao(À) 


Bulk modulus Bo (GPa) 


modulus Bo according to WC- 


GGA for Mgi-,.Mn,Se alloys Compound This work 
after structural optimization 
MegSe 5.94 
Mgo75Mno25$e 5.90 
M£go.50Mno.50$e 5.87 
Mgo25Mno.75$e 5.86 
MnSe 5.85 


Cal. Exp This work Cal. 

5.93 [38] 5.89 [32] 47.76 48.13 [38] 
5.92 [55] 49.00 [54] 
5.90 [36] 45.86 [36] 
- _ 49.72 - 

= = 52.17 - 

- — 55.60 - 

5.84 [56] 5.90 [54] 57.86 42.79 [56] 

130.88 [57] 


optimization by using the experimental values for binary 
compounds as of entered parameters. 

According to WC-GGA approximation, the calculation 
of these parameters for binary and ternary compounds 
in comparison to the experimental and other available 
theoretical results are presented in Table 1, where the lattice 
constants of the binary MgSe and MnSe are very close to 
the experimental values [32, 54] and stay in good agreement 
with theoretical calculations [36, 38, 55, 56]. 

Besides, there is some decrease in the lattice constants 
with increasing Mn concentration for Mg1_;Mn,Se (x — 
0.25, 0.5 and 0.75) because the ionic radius of Mn is smaller 
compared with Mg atom. While, the compressibility in this 
alloy increases when the concentration of Mn increase. 

We note that there is no experimental and theoretical cal- 
culation of lattice constants for Mg1_,Mn,Se compounds 
so our results constitute a prediction. As well, the computed 
lattice constants are used in next calculations of electronic 
and magnetic properties. 


3.2 Electronic Band Structures and Density of States 


Most of the physical properties of solids are related to 
the electronic band structures and the density of states, 
so the knowledge of these properties is important for 
its efficient use in optoelectronic, magneto-optic and 
electromagnetic devices. Here, by using the predicated 
lattice constant the spin-polarized band structures is 
computed for Mg1-,Mn,Se (x = 0.25, 0.50, 0.75 and 1) 
semiconductor with spin up (1) and down (|) orientations 
along high symmetry of the first BZ in the ferromagnetic 
zinc blende phase. The curves calculated within TB-mBJ 
approach to solve the problem of underestimation of band 
gap are displayed in Figs. 1, 2, 3, 4, and 5. Therefore in the 
case of x = 0 concentration we have treated and plotted the 
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band structure of the binary compound (MgSe) in the non- 
magnetic (NM) phase (Fig. 1) The Fermi level is set at the 
zero energy to be considered as a reference point. 

From the figures, for both spin orientation (up and down) 
with x = 0.25 to 1 the band structures reveals a direct band 
gap, where the minimum of the conduction band (CBM) 
and the maximum of the valence band (VBM) situated at l 
point of the BZ, which possess the semiconducting nature 
in these alloys. 

The direct band gaps of all alloys determined from the 
spin up case shown in Table 2 are plotted in Fig. 6. For 


Energy (eV) 


X WK 


W L F 


Fig. 1 The band structure of the binary MgSe (x = 0) in the non- 
magnetic phase 
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Fig.2 Spin-polarized band 
structures of majority spin (up) 


and minority spin (dn) for 14 = 
Mgo.75Mn025$e. The Fermi 12 K 7 
level is set to zero (horizontal 10 NZ 
dotted line) = DA 
8 ES 
$ DR 
: FE 
œ 2 
& 0 er 
d 2h — 
E Em | 
ui el 
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-6 
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further evaluation the band gap energy of Mg1_,Mn,Se can 
be depicted as a function of the manganese concentration 
and expressed by the following formula: 


Ec(Me,_.Mnsse) (x) = x Eg(Mnse) + (1 — X) Eg(Mgse) 


—x(1—-x)b (1) 


where E,(mg,_,Mn,se) (x) is the band gap energy of the 
Mgi-;Mn,Se, E,(Mese) is the band gap energy of the 


Fig.3 Spin-polarized band 
structures of majority spin (up) 
and minority spin (dn) for 
M£o.50Mn6.50$e. The Fermi 
level is set to zero (horizontal 
dotted line) 


Energy (eV) 
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CZ 


PANNE: 


MgSe, Ee(Mnse) is the band gap energy of MnSe, and b is 
the band gap bowing parameter of Mg1_,Mn,Se alloys. 
This equation gives b = 0.986 eV, according to the TB- 
mBJ approximation where the results shown in Fig. 6 was 
fitted by the expression (1) and follow the quadratic (2) 


Fr 


2 
ETB-mBn = 4094 — 3.281x + 0.986x 


E (2) 


Es 
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Fig.4 Spin-polarized band 
structures of majority spin (up) 
and minority spin (dn) for 
Mgo25Mn075$Se. The Fermi 
level is set to zero (horizontal 
dotted line) 


Energy (eV) 


As may be seen the energy gap decreases with increased 
concentration in our compounds. However, in the spin down 
channels, a symmetric band gap exists around the Fermi 
level (EF) for Mn-doped MgSe at all concentrations. 

In addition, the different band states of atomic and orbital 
origins can be explained by calculating the densities of 
states (DOS) and are displayed in Figs. 7, 8, 9, 10, and 11 
for x = 0, 0.25, 0.50, 0.75 and 1; the Fermi level EF is 
set to Zero and it is indicted by the vertical blue dashed 
line. 


Fig.5 Spin-polarized band TB-mBJ 
structures of majority spin (up) 14 NZ 
and minority spin (dn) for MnSe P d 
(x = 1). The Fermi level is set to 12 N 
zero (horizontal dotted line) 10 
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According to the crystal field theory, the PDOS curves 
show that the fivefold degenerate of 3d (Mn) states are 
divides in two parts: the threefold degenerate high-lying #4 
(dxy, dyz and dx) and the twofold generate low-lying eg (d,2 
and d,2_d,2) symmetry states, this is caused by the effect 
of the tetrahedral crystal field formed by surrounding (S, Se 
and Te) ligands; where the energies of e, states are lower 
than fn, states due to less Coulomb interaction. 

In the case of non-magnetic MgSe (x = 0) (Fig. 7), the 
DOS in the first range from — 11.5 to — 10.5 (eV) is mainly 
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Table 2 The calculated Es 


(spin-up band gap), total Band gap (eV) Magnetic moment (118) 
moments MT (8) per Mn, 
and the local magnetic Compound E ; TOP) Exp. Other work  MTt MM MMes M$ 
moments MM2/Mn/$e for each Our work 
concentration for 
Megi-;Mn,Se alloys (x = 0, MgsSe 412(NM)  4.00[31]  4.10[36] = _ = _ 
0.25, 0.50, 0.75 and 1) by using 
TB-mBJ functional +20 F8 
2.60 [36] 
Mgo.75Mn025Se 3.290 _ _ 5.000 4.383 0.007 0.027 
Mgo.50Mn0.50Se 2.682 _ _ 5.000 4.402 0.015 0.057 
Mgo25Mn075Se 2.257 _ _ 5.000 4.424 0.024 0.075 
MnSe 1.767 2.90 [58] 3.65 [57] 5.000 4.489 — 0.100 


dominated by ‘“4s” Se orbital. The following range from — 3 
(eV) to near the Fermi level (EF) contain “4p” orbital of Se, 
while the last region from 4 (eV) to 12 (eV) is dominated by 
“2p” “3s” Mg and 4p Se orbitals respectively. 

For x = 0.25, 0.50 and 0.75 (Figs. 8, 9 and 10) we 
found a great similarity behavior of the density electronic 
of states. While the lower part range between — 13.5 (eV) 
and — 11 (eV) is formed by 4sorbitals of Se atom for both 
spin up-down orientation. As well the DOS in the middle 
part range around — 5.5 (eV) to near the Fermi level EF 
comes essentially from the hybridization between the “3d- 
eg” “3d-t$" Mn and 4p Se, 3s Mg states in the spin up case, 
conversely the spin down case is dominated by 4p and 4s Se 
sates. The last one is located between 1.2 (eV) and 15 (eV) 
and divided into two sub-bands, the first one states from 1.2 
(eV) to 5 (eV) formed by 3d-e, and 3d-#2, Mn states for the 
spin down case. Where the second one is situated between 5 


Band-gap energy (eV) 


0,00 


0,25 0,50 


x composition 


0,75 


Fig.6 Band gap as a function of Mn concentrations for Mg,_,Mn,Se 
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(eV) and 15 (eV) mainly composed of mixture of 4s, 4p Se 
and 2p, and 3s Mg orbitals for both spin up and down. 

In the case of ferromagnetic MnSe (x = 1) (Fig. 11), the 
DOS in the first range from — 13.6 to — 12.7 (eV) mainly 
dominated by 4s Se orbital for both spin orientation (up 


TB-mBJ 
8 m 


! =— Total DOS 
! —— Total Mg 
— Total Se 


DOS(states/eV) 


DOS(states/eV) 


Energy (eV) 


Fig.7 The total and partial DOS of (2p, 3s) of Mg and (4s, 4p) Se for 
the binary MgSe (x = 0) in the non-magnetic phase. The Fermi level 
is set to zero (vertical dotted line) 
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TB-mBJ 
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È 0.25 
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; — Total Mn 
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Ë - Total Se 


DOS(States/eV) 


— 3d (Mn) 
— 3d-eg (Mn) 
: —— 3d-t2g (Mn) 
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Fig.8 Spin-polarized total and partial DOS of (2p, 3s) of Mg, (4s, 4p) 
of Se, and (3d, 3d-f2,, 3deg) of Mn in supercell of Mgo 75Mno25$e. 
The Fermi level is set to zero (vertical dotted line) 


and down). The second range from — 5.4 (eV) to near the 
Fermi level EF the DOS in the spin up case comes from the 
hybridization between 3d-e,, 3d-4 Mn and 4p Se orbitals, 
where the spin down case is dominate essentially by the 4p 
Se orbital. Also, the range between 0.5 (eV) and 5 (eV) is 
dominated by 3d-e, and 3d-#, Mn states for the spin down 
case, while from 5 to 15 (eV) the region is formed by 4p and 
4s Mg orbitals for both spin updown. 

On the other side for each concentration, no d-states 
peak of Mn at the Fermi level for both spin channels 
(up and down) which exhibit ferromagnetic semiconductor 
behaviors in these compounds is observed. 


3.3 Magnetic Properties and Exchange Coupling 


In this section, the values of total and local magnetic 
moment for DMSs Mgi-,Mn,Se at different concentra- 
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TB-mBJ 
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i---- 2p (Mg) 
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Energy (eV) 

Fig.9 Spin-polarized total and partial DOS of (2p, 3s) of Me, (4s, 4p) 


of Se, and (3d, 3d-#24, 3deg) of Mn in supercell of Mgo.59Mno.50$e. 
The Fermi level is set to zero (vertical dotted line) 


tions (x = 0.25, 0.5, 0.75 and 1) using TB-mBJ approx- 
imation were exerted by first principles calculation are 
summarized in Table 2. 

However, the magnetism in our selected compounds due 
to 3d (Mn) impurity where the Mn atom contributes two 
electronic to bonding of host valence band semiconductor 
resulting in Mn+? ion. As consequence the electronic 
configuration of 3d (Mn) in our alloys is = est3e) 
according to the Hund”’s rule for each concentration. Where 
all the spin-majority states are occupied but spin-minority 
states are either empty or partially occupied. 

Through Table 2, we observed that local magnetic 
moment of Mn doped MeSe at x = 0.25, 0.50, 0.75 and 1 is 
reduced less than the predicated Hund’s rule values Sug for 
total magnetic moment due to the p-d exchange interaction. 
On the other hand, a small contribution of local magnetic 
moment is induced on the nonmagnetic Mg and Se sites. 
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Fig. 10 Spin-polarized total and partial DOS of (2p, 3s) of Meg, (4s, 


4p) of Se, and (3d, 3d-f25, 3de,) of Mn in supercell of Mgo.25Mn0.75Se. 
The Fermi level is set to zero (vertical dotted line) 


After that, according to the mean field theory, the 
Hamiltonian given by the following equations [60, 61]: 


H = —Nop.S (3) 


(where No shows cation content, B expresses p-d exchange, 
while s and S, respectively, represent free hole and the Mn 
impurity spins) is employed to find tow important type of 
exchange constants Noa and NB which demonstrate the 
strength of the s-d and p-d couplings respectively. 

As well, the parameters can be calculated from the band 
structure and the magnetic properties by assuming the usual 
Kondo interactions [62]: 


Noa = RE (4) 
x(S) 

Mas (5) 
x(S) 
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Fig. 11 Spin-polarized total and partial DOS (4s, 4p) of Se, and (3d, 
3d-f28, 3deg) of Mn for the binary MnSe (x — 1). The Fermi level is 
set to zero (vertical dotted line) 


Here, AËc and AEYy are band-edge splitting of CBM and 
VBM respectively, calculated from given equations: 


(6) 
(7) 


_ pŸ d 
AËC = ECpmin — ECBMin 


AEv = Eng — AA 
x is the concentration of Mn atom and (S) in one 
half of magnetization per Mn ion. The predicted values 
of AECAEy, No and NoB respectively for each 
concentration of x are collected in Table 3. 

The negative sign of both N5B and Noa confirms the 
double-exchange mechanism in our compounds where the 
sd and p-d interactions are parallel and support FM 
character. In addition, we find that N5B is more negative 
than Noa this indicated that the exchange energy involved 
through the spin-down channel. 
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Table 3 Calculated 
conduction A Ec and valence Compound AËc AËy No NoB 
AEy band-edge spin-splitting 
and exchange constants Now, Mgo75Mn025$Se — 0.488 — 1.030 — 0.781 — 1.648 
NoB, for Mg1-,Mn,Se at x — Mgo.50Mn0.50$e — 0.791 — 1.771 — 0.632 — 1.417 
0.25, 0.50, 0.75 and 1 using 
k ; M M — 0.82 — 2. — 0.442 — 1.258 
TB-mBJ functional 80.25 Mno5$e me 0 ; 
MnSe — 0.228 — 2.986 — 0.091 — 1.194 
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In this work, we use the full potential linearized augmented plane wave (FP- 
LAPW#+lo) method under the framework of spin polarization density func- 
tional theory (SP-DFT) to study the structural properties of diluted magnetic 
semiconductors alloys, Mg:_,Cr,Se (x=0.25, 0.50 and 0.75), in the rock-salt 
ferromagnetic phase by using the WC-Cohen GGA approximation. To compute 
the electronic and magnetic properties of our compounds, we have used the 
Tran Blaha modified Becke-Johnson potential (TB-mBJ) combined with PBE 
approximation. The results of electronic properties show that our compounds 
have a half-metallic behavior with a spin polarization of 100% at the Fermi 
level. Moreover, at all concentrations, the magnetic moment has been esti- 
mated as equal to 4.0 (18). Furthermore, to validate the effects resulting from 
the exchange splitting process, we calculate the values of the spin-exchange 
constants Noa and NB, respectively. Finally, we present in detail the effects of 
the temperature and the pressure on some macroscopic thermodynamic 
parameters by using the quasi-harmonic Debye model. 


Key words: Ab initio calculation, electronic structure, magnetic properties, 
spintronic, TB-mBJ, DMS compounds, thermodynamic 


properties 


INTRODUCTION 


There has been huge progress in the area of semi- 
conductors, and many works are devoted to designing 
new materials; today the search continues to improve 
miniaturization, as well as the speed ofthe components, 
which highlights the limits of classical electronics. That 
is why more projects are studying a new type of elec- 
tronics: spintronics or spin electronics, which use the 
spin of the electron in addition to its charge to guide the 
electrons, store information and create currents. 

Spintronics continues to develop in many direc- 
tions, and researchers search for new materials to 
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meet the technology needs in various fields of 
applications. In this context, diluted magnetic 
semiconductors (DMS) or semimagnetic semicon- 
ductors have attracted special interest of both sci- 
entists and technologists during the last few years 
because they opened a way for the development of 
new physical properties on bulk and two-dimen- 
sional (2D) materials. 

In our case, these materials can be obtained by 
doping a classical semiconductor (II-VI, IV-VI, III 
V, IV, and II-VI-VL) with magnetic elements such 
as transition metals (TM) or rare earth (RE),! which 
introduces freedom of spin into semiconductor de- 
vices. 

Moreover, DMS compounds are good candidates 
to obtain half-metallic (HM) appearance, as they 


128 


A Systematic Study of Structural, Magneto-Electronic and Thermodynamic Properties of Mg:_,Cr,Se 3249 
DMS Alloys in the Rock-Salt Phase for Spintronic Applications 


have metallic characteristics for one spin direction 
and are semiconductors for the other spin direction, 
while at the same time exhibiting a polarization of 
100% at the Fermi level (Er).? 

We recall that, nowadays, HM in DMS com- 
pounds is a key ingredient for the future high-per- 
formance spintronic devices, such as light-emitting 
diodes, magnetic sensors, ultra-fast optical 
switches,” spin-polarized transistors,‘ spin valves 
and magnetic read head sensors of computers.” In 
addition, the thermal properties of crystalline 
materials are very important and useful in many 
industrial applications involving high pressure and 
high temperature. In this direction, numerous re- 
search groups depicts the thermodynamic proper- 
ties of half-metallic compounds” to show their 
multifunctional characteristics. 

In fact, other theoretical results!°-!” have pre- 
sented the impact of the doping of binary semicon- 
ductors with magnetic transition elements in 
several crystalline phases with a large variety of 
concentration to predict their usage in spintronic 
applications. 

On the other hand, magnesium chalcogenides are 
technically remarkable elements of II-VI semicon- 
ductors, having a wide range of applications ranging 
from catalysis to microelectronics, * and particu- 
larly they are a good choice for protective coatings 
due to their hardness, high melting point, high 
thermal conductivity, and large bulk modulus. !? 

As a class of II-VI semiconducting materials, 
magnesium selenide (MgSe) is a very attractive 
binary semiconductor with both n- and p-type dop- 
ing ability/?1 and could be used in high-tempera- 
ture, high-power blue and ultraviolet wavelength 
optics”? ?{ in ZnSe-based laser diodes (LDs) and 
light emitting diodes (LEDs) as a cladding mate- 
rial'$ and in spintronic devices.?°?° 

However, the crystallographic structure of bulk 
MgSe exists in different forms: zinc blende (B3), 
wurtzite (B4), and rock-salt (B1).!° Each structure 
is different from another on the basis of lattice 
parameter and structure constants. 

The structure of zinc blende is a face-centered 
cubic (fcc) lattice with diatomic basis, while that of 
wurtzite is symmetrically hexagonal with four 
atomic bases. Meanwhile, the MgSe rock-salt 
structure is cubic and can be found in the 
nanocrystalline form! and has an excellent 
technological application.®? This is the reason be- 
hind our study of the rock-salt structure in detail as 
most of the interesting properties of the bulk 
materials are related to nanostructures. 

In addition, doping rock-salt structures with 
certain impurities may change and enhance their 
properties,” and with this background, it is 
favorable to study our semiconductor with magnetic 
doping effects to explore the possibility of HM 
behavior in this structure. 

The objective of the present work is to investigate 
the structural, electronic, magnetic and thermody- 


129 


namic properties of Mg:_,Cr,Se (x=0.25, 0.50, and 
0.75) DMS compounds in the ferromagnetic rock- 
salt phase using ab initio calculations and Gibbs 
quasi-harmonic Debye model. 

Our paper is organized as follows: we shall de- 
scribe our computational detail in the next section 
and the third section is dedicated to results and 
discussion. Finally, the summary is given in the 
fourth section. 


COMPUTATION METHOD 


Our calculations are performed with spin-polariza- 
tion density functional theory (SP-DET)” using 
WIEN2k*% package under the framework of the full 
potential linearized augmented plane waver (FP-LAPW 
+10)? method to solve the Kohen-Sham equation.*® 

However, to simulate the properties of MgSe 
semiconductor doped with Cr impurity, we have 
used a (1x 1x1) supercell with eight atoms of the 
rock-salt (RS) structure Fm3m (SG: 225) with a 
cubic symmetry Fig. 1, four atoms are of Mg and 
four of Se. The structure of Mg:_,Cr,Se is obtained 
by replacing the Mg atoms with Cr to get the com- 
pounds with x=0.25, 0.50, and 0.75 concentrations, 
respectively. 

In this case, the generalized gradient approxi- 
mation function developed by Wu and Cohen (WC- 
GGA)*° is taken as an exchange-correlation (XC) of 
the electrons to represent the equilibrium struc- 
tural parameters. As it is well known, the Tran- 
Blaha modified Becke Johnson (TB-mBJ)“ ap- 
proach is used for the precise calculations of the 
electronic structure, band gap semiconductor and 
materials based transition metal. #14 

Hence, we propose to use the TB-mBJ approach 
combined with Perdew-Burk-Ernzerhof approxima- 
tion‘ in the form of TB-mBJ_PBE to predict the 
electronic and magnetic properties of Mg:_,Cr,Se 
compounds. 

In our case, the plane wave expansion with R,,+x 
Kax Was chosen to be equal to 8, where R,,4 is the 
smallest radiïi of muffin-tin spheres and K,,.. is the 
cut-off for the wave function basis. The Mg: 35°, Cr: 
3d°4s! and Se: 3d!°4s?4p* electrons are treated as 
valence electrons. The atomic radius (muffin-tin 
RMT) was chosen equal to 2.31, 2.20 and 2.31 
atomic units (a.u) for Mg, Se, and Cr, respectively, 
whereas the Fourier expanded charge density was 
truncated at Gnax=12 (a.u) !, and the k-mesh of 12 
x 12 x 12 equivalent to 2000 k points was selected 
for the Brillouin zone (B2). For each calculation, the 
total energy convergence is equal to 107* (Ry) when 
the charge of the system is converged to 0.0001e. 

Finally, the results obtained of the ground state 
parameters from the Æ (V) curves with WC-GGA 
approximation in the ferromagnetic phase, are used 
to investigate the thermodynamic properties of Mg. 
Cr,Se (x=0.25, 0.50, and 0.75) compounds through 
the quasi-harmonic Debye model as implemented in 
the Gibbs program.“° 


Bordjiba, Meddour, and Bourouis 


Fig. 1. The (a-d) crystal structures of Mg:_,Cr,Se (x=0.0, 0.25, 0.50, and 0.75). The gray, blue, and light-yellow spheres, respectively, indicate 
Cr, Mg, and Se atoms (Color figure online). 


Table I. Calculated lattice constant ao, bulk modulus B,5, formation energy E;, and cohesive energy Ec 
according to WC-GGA for Mg,_,Cr,Se compounds x=0, 0.25, 0.50 and 0.75 in the rock-salt phase 


Lattice constant & (À) Bulk modulus B, (GPa) E; (eV/atom) E, (eV/atom) 
Compound This work Cal. Exp This work Cal. Exp This work This work 
MgSe 5.46 5.46% 5.47% 64.62 65.0'° : : _- 
5.51 61.0 
Mgo.75Cro.25Se 5.41 — = 67.69 = = — 0.657 — 2.928 
Mgo.50Cro.509e 5.36 cu = 73.64 RÉ = — 0.528 — 3.400 
Mgo.25Cro.75Se 5.31 LL Eu 76.46 = — — 0.398 — 3.870 


RESULTS AND DISCUSSION 
Structural Properties 


In order to obtain the ground state properties of 
Mg:_,.Cr,Se compounds in the ferromagnetic rock— 
salt structure for each concentration (x=0, 0.25, 
0.50, and 0.75), we apply the WC-GGA*° approxi- 
mation to find the structural properties such as the 
equilibrium lattice constant (ao) and the bulk mod- 
ulus (Bo), where the total energy is fitted with 
Birch-Murnaghan equation of states.“ The calcu- 
lated structural parameters are displayed in 
Table I. 

Thus, the value of the lattice parameter for the 
binary MgSe is in good agreement with other theo- 
retical works!” and equal to the experiment va- 
lue.*° In addition, the composition dependence of 
lattice constants as a function of Cr-doping concen- 
tration for the studied compounds shows that (@o) 
decrease with the increase of x, this arises because 
the atom of Cr is smaller than Mg. 

On the other hand, the compressibility in these 
alloys increases when the concentration of Cr in- 
crease. We recall that there is no experimental and 
theoretical calculation for the structural properties 
of Mg:_,Cr,Se (x=0.25, 0.50, and 0.75) compounds; 
therefore, we have used the computed lattice con- 
stants in the next calculations. 

In addition, to justify the existence of our com- 
pounds Mg:CrSez, Mg2Cr2Se4, and MgCr;$e, in the 
rock-salt ferromagnetic structure, we have com- 


puted the defect formation energy Et. It is usually 
used to measure the phase stability in solid 
states“ °1 by using the following equation: 


Ef = [Erota(Mg,CriSec) — aEbulk (Mg) — bEbux(Cr) 
—CErux(Se)] 
(1) 


where ErtAa(Mg.Se,Cr;) represents the total energy 
of the compound state, Ekun(Mg), Ebunx(Cr), and 
Epux(Se) are the total energies per atom of fully bulk 
Mg, Cr, and Se respectively, and a, b, c are the num- 
bers of Mg, Cr, and Se atoms in the unit cell formation. 

Furthermore, based on the previous calculation 
method in other systems," we have also 
premeditated the cohesive energy Æ, to assess the 
thermodynamic permanence of our compounds 
Mg:CrSes, Mg°CrSe4, and MgCr3Se4 by employing 
the following relation: 


Ec = (Erota( Mg, CriSec) — GE atom(M£) 


—DE atom(Cr) — CEatom(Se)] Le 


Here E tom is the energy of the free atom in its 
ground state and calculated from each isolated atom 
in a cubic cell with a 15 (A) lattice constant to ne- 
glect the interactions between neighbor atoms. 

We realized that both defect formation and cohe- 
sive energies show negative values as given in 
Table I. These results indicate the presence of our 
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compounds in the rock-salt structure, which sug- 
gests that the system under study can exist in this 
phase. 


Electronic and Magnetic Properties 


Spin Polarization Band Structure and Density 
of States 


In this section, we study the electronic properties 
of MgSe doped with 3d transition metal chromium 
(Cr) at x=0.25, 0.50 and 0.75 concentrations in the 
ferromagnetic rock-salt phase. The knowledge of 
these properties in the DMS compounds provides 
valuable information for their potential utility in 
spin electronic devices. 

By using the theoretical lattice constants and the 
TB-mBJ_PBE approach, we have computed the to- 
tal density of states (TDOS) and the band structures 
for both spin channels in up (f) and down (|) ori- 
entations of Mg,_,Cr,Se along high symmetry of the 
first Brillouin Zone (B2). All the results are pre- 
sented in Figs. 2, 3, and 4, where the Fermi level 
(Er) is set at zero energy to be considered as a ref- 
erence point. 

For these materials, both the band structure and 
the total density of state of the majority spin (up) 
show a metallic nature due to the presence of a large 
peak centered at the Fermi level. On the other 
hand, the minority spin (down) accentuates a direct 
band gap semiconductor with the top of the valence 
band (TVB) and the bottom of the conduction band 
(BCB) located at T point of the BZ, giving a gap of 


2.667 (eV), 2.044 (eV) and 1.831 (eV) for 
Mgo75Cro25Se, Mgo.50Cro.50Se and Mgo25Cro.75Se, 
respectively. 

It is also evident from Figs. 2, 3, and 4 that the 
band gap can also be decreased by increasing the 
concentration of Cr in alloys. In addition, the Fermi 
level does not lie in the middle of the impurity in the 
total density on the spin down case, suggesting that 
these three DMS compounds exhibit half-metallic 
character with 100% spin polarization at Fermi le- 
vel. 

Moreover, The HM gap Gym is a very significant 
parameter to determine the importance of DMS”° in 
the spintronic devices; it is defined as the minimum 
of the lowest energy of majority-spin and minority- 
spin conduction bands with respect to the Fermi 
level, and the absolute values of the highest energy 
of majority-spin and minority-spin valence bands. 
The calculated values of HM gap are listed in 
Table II, where it can be seen that the HM gap 
values become smaller with increasing Cr concen- 
tration up to x=0.75. 

In order to study the atomic contribution to the 
ferromagnetism in more detail, the partial density 
of states (PDOS) for the Mgz5CrossSe, 
Mgo5oCro5oSe, and Mgo25Cro.75Se compounds have 
been plotted in Figs. 5, 6, and 7. The Fermi level Er 
is set to zero and indicated by the vertical red da- 
shed line. 

The spin polarization of a system at (Ex) is de- 
fined as the following expression: 


TB-mBJ_PBE 


M9 75C 0255 


SEVRES. 


Energy (eV) 


Half-Metallic! 


a | : L AN 
—_ SE J 
Semi-conductor 


E,=2.667 eV 


3 6 
Density of States (States/eV) 


Fig. 2. Spin-polarized band structures and the total density of states of majority spin (up) and minority spin (dn) for Mgo 75Cro 25Se by using TB- 


mBJ_PBE. The Fermi level is set to zero (horizontal dotted line). 
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Fig. 3. Spin-polarized band structures and the total density of states of majority spin (up) and minority spin (dn) for Mgo 50Cro.50Se by using TB- 
mBJ_PBE. The Fermi level is set to zero (horizontal dotted line). 
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Fig. 4. Spin-polarized band structures and the total density of states of majority spin (up) and minority spin (dn) for Mgo25Cro.75Se. The Fermi 
level is set to zero (horizontal dotted line). 


N\(Er) - NU(Er) where N'(Er) and N!{Er) are the spin up and down 
(3) densities of states at Er. 
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Table II. Calculated total moments M°°* (up) per Cr, and the local magnetic moments MVC" for each 
concentration for Mg,_,Cr,Se compounds (x=0.25, 0.50, 0.75) by using TB-mBJ_PBE functional in the rock- 


salt phase 


Compound E, ‘P°Y") (eV) Guam (eV) MT (up) M (8) ME (uB) M (8) 

Mgo.75 Cro.25Se 2.667 0.902 4.000 3.654 0.006 —0.026 

M£go.50 Cro.soSe 2.044 0.382 4.000 3.639 0.013 —0.032 
TB-mBJ_PBE TB-mBJ_PBE 
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Fig. 5. Spin-polarized total and partial DOS of (2p, 3s) of Mg, (4s, 
4p) of Se, and (34, 3d-4,, 34.) of Cr in a supercell of 


Mgpo.75Cro.25Se. The Fermi level is set to zero (vertical dotted line) by 
using TB-mBJ_PBE. 


In Figs. 5, 6, and 7, we observe that for all three 
compounds the lower part of the valence band (VB), 
between —15 (eV) and —-12 (eV) is occupied by 4s (Se) 
states for both spin cases and the DOS in the middle 
part of VB, in the energy range from —7 (eV) to -8.5 
(eV), is formed by 4p (Se) and 3d (Cr) states. 
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Fig. 6. Spin-polarized total and partial DOS of (2p, 3s) of Mg, (45, 
4p) of Se, and (34, 3d-t,, 34.4) of Cr in a supercell of 
Mgo.50Cro.508Se. The Fermi level is set to zero (vertical dotted line) by 
using TB-mBJ_PBE. 


In addition, according to the crystal field theory, 
the octahedral crystal 3d Cr states split in fivefold 
degenerate, low-lying {2 (dxy,dyz, and dxz) and 
high-lying eg (d2 and d,2 — dy) symmetry states. 

Furthermore, the values of the density of state at 
Er in the spin up case is equal to 
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Fig. 7. Spin-polarized total and partial DOS of (2p, 35) of Mg, (45, 
4p) of Se, and (34, 3d23, 3dg) of Cr in a supercell of 
Mgo25Cro.75Se. The Fermi level is set to zero (vertical dotted line) by 
using TB-mBJ_PBE. 


N'(Er) = 0.61 (states/eV), N'(Er) — 1.13 (states/eV) 
and N'(Er) = 1.86(states/eV) per Cr states for 
Mgo.75Cro.25Se, Mgo.50Cro.soSe, and Mgo:25Cro.75Se 
receptivity, where most of this density comes from 
the Cr-3d-e, states. We note that the contributions 
from Se-4p and Cr-3d-t2, states are small. We recall 
that the metallic character in the spin up case of 
these compounds results from the large hybridiza- 
tion between 4p (Se) and 3d (Cr) states, and we re- 
call again the density of states at Er is equal to zero 
in the spin down case (N{Er#)=0) at all concentra- 
tions, which gives a polarization P=1. 

In addition, the conduction band (CB) ranging 
from 6 (eV) to 12.5 (eV) is principally formed by 3s 
and 2p Mg states in both spin channels. 

We have also calculated the exchange splitting 
energy (A(d)=d,-d;), which is defined as the sep- 
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aration between the spin up and down peaks, and 
we compared it with the anionic crystal filed energy 
(Ecrys=dr2 deg) to explore the existence of a 
stable ferromagnetic state. Also, the p-d energy 
interactions (A,({pd)=E}-E?) of the valence bands 
top for the spin up and down orientations was pre- 
sented to understand more about ferromagnetism in 
our compounds. 

According to Table III, the calculated values of d— 
d coupling are matched higher than the E,,.,.=d;29- 
deg. These results confirm the presence of ferro- 
magnetism in these compounds, and the negative 
exchange-splitting gap A,(pd) values mean that the 
effective potential for the minority spin is more 
attractive as compared to the majority spin, which 
characterizes the spin-polarized systems. 

From these results, we confirm the HM behavior 
in our compounds, which can be applied as elec- 
trodes”® and exploited in many spintronic tech- 
nologies such as giant magnetoresistance (GMR), 
tunnel magnetoresistance (TMR), random access 
magnetic memory (MRAM), and quantum qu- 
bits," 


Electron Charge Density 


In order to study the bonding properties between 
atoms, we have examined the electron charge den- 
sity by implementing TB-mBJ_PBE approximation 
for the considered structures in the (110) plane by 
talking into account the spin up and spin down 
configurations along with 2D bonding charge den- 
sity distributions, as shown in Fig. 8. It is observed 
from the figure that there is apparently no notice- 
able difference in the total valence charge densities 
at all concentrations, and the contour plot shows the 
ionic character of Mg-Se bond due to the Pauling 
electronegativity difference and a strong covalent 
Cr-Se bond, while the charge distribution among 
the atoms indicates the presence of partial covalent 
character. In other words, the valence charge 
around the Mg sites is due to s states, and both s 
and p states are responsible for the charge around 
the Se atoms. 


Magnetic Moment and Exchange Coupling 


By using the TB-mBJ_ PBE approximation, we 
have presented the total and the local magnetic 
moment of Mg:.,.Cr,Se DMS compounds in the 
rock-salt ferromagnetic phase for x=0.25, 0.50 and 
0.75. The obtained results are summarized in 
Table IT. 

Further, according to the hand’s rule, the elec- 
tronic configuration of 3d (Cr) in Mg;_,Cr,Se for 
each concentration is (de), hence the majority 
spin of 8d (Cr) is partially filled by three electrons 
whereas the minority spin is empty. 

From the Table II, the local magnetic moments of 
Cr impurity such as 3.654, 8.639 and 3.668 (18) at x 
=0.25, 0.50 and 0.75 concentrations respectively are 
smaller than the predicted Hund’s rule values of 4 
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Table III. Calculated exchange splitting A,(d), crystal filed energy E..,, band-edge spin-splitting A,(pd), and 
exchange constants Noo, Nof for Mg:_,Cr,Se compounds at x=0.25, 0.50 and 0.75 


Compound AE crystal A, (d) À; (pd) Not Nob 

Mgo7s Cro25Se 0.463 83.918 —0.514 8.917 — 1.028 
Mgv.50 Cro.5oSe 0.490 4.000 —0.706 1.810 —0.706 
Mgo25 Cro7sSe 0.816 4.544 — 0.680 1.012 —0.453 


1,000 
0,9600 
0,9200 
0,8800 
0,8400 
0,8000 
0,7600 
0,7200 
0,6800 
0,6400 
0,6000 
0,5600 
0,5200 
0,4800 
0,4400 
0,4000 
0,3600 
0,3200 
0,2800 
0,2400 
0,2000 
0,1600 
0,1200 
0,08000 
0,04000 
0,000 


TB-mBJ_PBE 


"y 


Mg, Cri 75e 


ZO 


Fig. 8. Electron charge density of (a) Mgoz5Cro25Se, (b) Mgo5oCro.5oSe, and (c) Mgo25Cro.75Se compounds. For both spin up and down 
configurations computed by using TB-mBJ_PBE. 
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(8) for total magnetic moment which they caused 
by the hybridization between d-states of Cr and p-Se 
states exchange interaction. 

In addition, we can see a small local magnetic 
moment localized on the nonmagnetic (Meg, Se) sites. 
It is noteworthy that the negative local magnetic 
moments sign of Se indicate the anti-ferromagnetic 
interaction between Cr magnetic spin of Se anion, 
while the positive values of magnetic moment at the 
Mg sites indicate the ferromagnetic interaction be- 
tween Mg and Cr spins. 

In addition, using the Hamiltonian equation gi- 
ven by the following formula®!: 


H = -Noip,.S (4) 


(where N; shows cation content, B, expresses p-d 
exchange, while s and S represent, respectively, free 
hole and the Cr impurity spins) we can calculate two 
important parameters: the s—-d exchange constant 
Noa (conduction band) and p-d exchange constant 
NoB (valence band) from the magnetic properties 
and the band structures by assuming the usual 


Kondo interactions.°152 
ÂAËc 
Not x(S) (5) 
_ AËEy 
Nof = xtS) (6) 


Here, AËc and AËy are the band-edge splitting of 
TVB and BCB respectively, calculated from given 
equations: 


AËEc = Ejve ou Elve (7) 


AEy — Eice =, Eice (8) 


x is the concentration of Cr atom and (S) in one half 
of magnetization per Cr ion. 

The predicted values of AEc, AEËv, Noû, and NB 
for each concentration of x are listed in Table III. It 
shows that values of Noa are positive, whereas Nof 
values are negative for all three compounds indi- 
cating that the anti-ferromagnetism and ferromag- 
netism occur between the valence band and 3d 
states of Cr and also between the conduction band 
and Cr valence spins, respectively. 


Thermodynamic Properties 


To get more information about the thermody- 
namic properties in the ferromagnetic rock-salt 
phase of Mg:_,Cr,Se (x=0.25, 0.50, and 0.75) DMS 
compounds in the range of 0-20 (GPa) pressures, we 
applied the quasi-harmonic Debye theory of crystal 
as implemented in Gibbs program“%%%° with the 
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FP-LPAW method to investigate the macroscopic 
parameters as a function of temperature in the 
range of 0—-1200 (K). 

Figures 9 and 10 display the dependence of heat 
capacity C, and Cy with temperature and pressure 
of Mg:_,Cr,Se. It can be seen that from 0 to 300 (K) 
Cy and CE increase exponentially following the T' 
law, and have similar behaviors for the three com- 
pounds. 

For T'>300(K), Cv tends towards a limit value 
known as the Petit-Dulong limit, but the behavior of 
Cp is different than Cy because it does not tend to a 
constant value at high temperature. 

The Cy values at room temperature and zero 
pressure of Mgoz5Cro25sSe, MgosoCro5oSe and 
Mgo5Cro25Se compounds are, respectively, 
69.0025 (J mol ! K !), 69.34634 (J mol ! K !) and 
72.1579 (J mol ! K !). 

The entropy S curves as a function of temperature 
at various pressures of Mg:_,Cr,Se are plotted in 
Fig. 11. The curves show that S increases sharply 
when the temperature increases and decreases with 
the increase of the pressure at a given temperature. 
As the temperature increases, the vibrational con- 
tribution to the entropy increases and, therefore, 
the entropy increases with temperature at different 
concentrations of the doping element (Cr). 

In Fig. 12. We present the relationship between 
the volume and the temperature at different pres- 
sures for our DMS compounds. The results show 
that at a given pressure, V increases slowly with the 
temperature, especially at high pressures. On the 
other hand, V decreases with the increase of both 
pressure and the concentration of the doping ele- 
ment (Cr). 

On the other side, the thermal expansion coeffi- 
cient a is an essential parameter for predicting the 
thermodynamic equation of state. Figure 13 depicts 
the effect of the temperature and pressure on this 
coefficient for the Mg:_,Cr,;Se compounds at differ- 
ent concentrations. For a given pressure, «à in- 
creases quickly with temperature up to 300 (K), and 
above this temperature it tends to a nearly constant 
value. As already seen on the curves of the Fig. 12, a 
decreases with the increases in both concentrations 
of the doping element (Cr) and the pressure. 

Figure 14 gives the evolution of the bulk modulus 
B of Mg:_,Cr,Se compounds with the temperature 
at different pressures. 

We observe that, at given pressure, B decreases with 
the temperature and when the values of pressures in- 
crease B increase. On the other hand, itis clear that the 
increase of the concentration in the compounds does 
not have a significant effect on the B values. 


CONCLUSION 


Diluted magnetic semiconductors (DMS) have 
both semiconducting and ferromagnetic properties, 
which provides the possibility to manipulate elec- 
tron charge and spin at the same time. In this pa- 
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Fig. 9. Variation of the heat capacities CR versus temperature at different pressures of (a) Mgoz5Cro25Se, (b) Mgo5oCro.5oSe, and (c) 
Mgo25Cro.75Se compounds by using WC-GGA. 
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Fig. 10. Variation of the heat capacities C, versus temperature at different pressures of (a) Mgoz5Cro2sSe, (b) Mgos5oCrosoSe, and (c) 
Mgo25Cro.75Se compounds by using WC-GGA. 


(a) (b) (c) 
Mg, Cr, Se] 
200 LM orsT025 + 200 200 
= 150} = 150 = 150 
ES | 5 5 
Ë Ë Ë 
2 E 14 3 2 
=: 100 0 GPa — 100 = 100 
è | —+—46GPa | 2 ES 
£ —+— 8 GPa © © 
£ —+—12GPa| € £ 
u L 4 [mr 
50 16 GPa W 50 50 
—— 20 GPa 1 —— 20 GPa 
0 0 ñ ni 1 1 ñ 1 … L ñ 0 
0 200 400 600 800 1000 1200 0 200 400 600 800 1000 1200 0 200 400 600 800 1000 1200 
Temperature (K) Temperature (K) Temperature (K) 


Fig. 11. Variation of the entropy S versus temperature at different pressures of (a) Mgo75Cro25Se, (b) Mgo.5oCro.5soSe, and (c) Mgo 25Cro 75Se 
compounds by using WC-GGA. 


per, we used the ab initio calculation and FP-LAPW From our study, we conclude that when the 
+lo method to investigate the structural, electronic, chromium concentration increases (x=0.25, 0.50 
magnetic and thermodynamic properties of Mg: and 0.75) the lattice parameter decrease. The band 
_.Cr,Se compounds in rock-salt phase. structures and the density of states show that the 
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Fig. 12. Variation of Volume V versus temperature at different pressures of (a) Mgoz5Cro25Se, (b) Mdo5oCro.5o0Se, and (c) Mgo25Cro.75Se 
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Fig. 13. Variation of the thermal expansion coefficient a versus temperature at different pressures of (a) Mgo 5Cro25Se, (b) Mdo 5oCto508e, and 
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Fig. 14. Variation of the bulk modulus B versus temperature at different pressures of (a) Mg xCro2sSe, (b) Mdo5oCto50Se, and (c) 


Mdo25Cro.75Se compounds by using WC-GGA. 


spin polarization in the rock-salt MgSe semicon- 
ductor can be introduced via the Cr doping and 
gives half-metallic behavior for each concentration, 
which this investigation predicate that our com- 


pounds are a potential candidate in the future 
industrial spintronic devices. 

For the three compounds, a large magnetic mo- 
ment is observed on the magnetic atomic sites with 
a small local magnetic moment on the non-magnetic 
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Mg and Se atomic sites caused by hybridization of p 
(Se) and d (Cr) states. 

The results provided by the quasi-harmonic De- 
bye-model, implemented in the Gibbs code, have 
shown that the concentration of the dopant (Cr) and 
the pressures does not have a significant impact on 
both C, and Cy. On the other hand, we observe an 
effect on S, V, and «. 

Finally, in the absence of the experimental data, 
this theoretical study can be considered as a good 
prediction and will be a reference for future work on 
these materials. 
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